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ABSTRACT ^ 

Was PSl-lOjh (jOezari et al.l 1201Z ). an optical/UV transient discovered by the Pan- 
STARRS Medium Deep Survey, the tidal disruption of a star by a massive black hole? 
We address two aspects of the problem: the composition of the putative disrupted 
object (using the spectroscopic data), and the energetics of the observed gas and 
radiation (using the photometric data). We perform photoionization calculations and 
compare with the observed lower limit of the line ratio LHeii 4686 /^Ha > 5 to argue 
that this event was not the disruption of a solar-ty pe star , and i nstea d was likely the 
disruption of a helium core (as first proposed by iGezari et al.1 l2012ll . Disruption of 
such a dense object requires a relatively small central BH, Mbh ^2x 1O 5 M 0 . We use 
the photometric data to infer that PSl-lOjh comprised an outflow of ~ O.OIMq of gas, 
escaping from the BH at ~ 1000kms -1 , and we propose that this outflow was driven 
primarily by radiation pressure trapped by Thomson and resonance line scattering. 
The large ratio of radiated energy to kinetic energy, E Ta d/EK ~ 10 4 , together with 
the large value of E la ^ ~ 2 x 10 51 erg, suggests that the outflow was shocked at large 
radius (perhaps similar to super-luminous supernovae or the internal shock model 
for gamma-ray bursts). We describe puzzles in the physics of PSl-lOjh, and discuss 
how this event may help us understand future tidal disruptions and super-Eddington 
accretion events as well. 

Key words: galaxies: nuclei - black hole physics 


1 INTRODUCTION 


On M ay 31 , 2010, the Pan-STARRS Medium Deep Survey 
llKaiser e t al .1120 1 ()l) first dete cted a remarkable event called 
PSl-lOjh ( Gezari et alJl20l3 l. Over the next two months, 
the event rose in brightness by a factor of ~ 200, then 
declined over the subsequent year. Follow-up observations 
showed colors and spectra inconsistent with a supernova 
(SN) or active galactic nucleus (AGN). These authors in¬ 
terpreted the event as the tidal disruption of a helium-rich 
stellar core. 

Tidal disruption events (TDEs) must happen every so 
often in the centers of galaxies, when a star passes so close to 
the central massive black hole (BH) that the BH’s tidal grav¬ 
ity exceeds the binding gravity holding the star together. 
The star disrupts, sending a fraction of stellar gas on el¬ 
liptical orbits that eventually return to the BH; the stream 
of stellar gas shocks on itself and begins to accrete onto 
the BH, releasing a flare of electromagnetic energy that we 


in pr i nciple can observe (e.g., lLacv. Townes fe Hollenbachl 
Il982l : [Reesll 19888 1. TDEs are interesting in their own right, 
and also as signposts of otherwise quiescent BHs, probes of 
stellar dynamics in the nuclei of galaxies, and informants 
about accretion physics. 

A few dozen transient events have been identified as 
candidate TDEs over the last c ouple of decades (e.g., 
iKomossal 120021 : 1 Gezari et al.l 120091 : iBloom et al.l 120111 1 , but 
the advent of wide-held transient surveys is now dramati¬ 
cally improving t he coverage of data available for any indi- 
vidual event (e.g.. IGezari et al .ll2oi~3 : IChornock et al.ll2014l : 

Ar cavi et al. 120141 1. PSl-lOjh is special because: 

(i) It is the first event whose light curve follows the rise 

to peak brightnesLI (not just the later-time decline), 

(ii) It is the first event whose optical spectrum shows 

lines at all (with the possible exception of TDE2; 
Ivan Velzen et ah IB), 

(iii) The spectrum shows He II emission but no hydrogen, 


* UBC Science Teaching and Learning Fellow. Email: 
linda@phas. ubc. ca 


1 PSl-lOjh is also only the second optica lly-selected event t o be 
followed up in real time, after PTFlOiya dCenko et al.|[2012h . 
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suggesting that the event was the disruption of the core of a 
red giant star rather than the disruption of a main sequence 
star, and 

(iv) As a consequence of (3), the host galaxy may harbor 
a particularly low-mass central BH. 

In addition, ultraviolet observations of PSl-lOjh indicate a 
relatively cool effective temperature T e a ~ 3 x 10 4 K that 
appears to have remained mysteriously constant for most of 
a year, while the luminosity ranged by a factor of ~ 50. 

We seek to develop a model that explains the observa¬ 
tions and offers insights into the tidal disruption process and 
accretion physics. 


1.1 Observed properties of PSl-lOjh 

Figure 2 in iGezari et al.l ((20121 ') shows the observed light 
curve of PSl-lOjh in optical (gr, r, i, z from Pan-STARRS) 
and ultraviolet ( near-UV from the Galaxy Evolution Ex¬ 
plorer: GALEX: iMartin et al.l 1200 -tI ). The event reached a 
maximum g-band luminosity of 2 x 10 43 ergs -1 on July 12, 
2010; after peak, the event was detectable to Pan-STARRS 
for a little over a year (until August 2011), over which time 
the gr-band luminosity faded by a factor of ~ 50. GALEX ob¬ 
served the event about 20 days before g-band peak and again 
about 250 days afte r pe ak. For th e two epochs that GALEX 
data was available, IGezari et alJ J2OI2I) fit the optical and 
NUV points of the spectral energy distribution (SED): they 
find that the SED at both epochs can be fit by a blackbody of 
temperature T = 2.9 x 10 4 K, assuming a galactic extinction 
of E(B — V) = 0.013 mag. This temperature would imply 
a peak bolometric luminosity of Lbol ~ 2.2 x 10 44 ergs -1 . 
The extinction may actually have been greater than this, 
leading to a hotter blackbody temperature; the limit of 
E(B — V) < 0.08 mag from the He II line ratio (discussed 
belo w) implies a t empe rature of 5.5 x 10 4 K. 

IGezari et alJ (120121 ) followed up the event with an opti¬ 
cal spectrum from MMT about 20 days before (17-band) peak 
and again about 250 days after peak. The optical spectra 
can be fit by stellar light from the host galaxy, plus a black¬ 
body continuum from the event (T = 2.9 x 10 4 K), and 
intriguingly—broad He II emission lines. In the spectrum at 
time tspect = ipeak — 22 days, two He II lines are visible, at 
rest wavelengths A = 4686 A (n = 4 —> 3) and A = 3203 A 
(n = 5 —> 3). The line at 4686 A has a linewidth of 
9000 km s -1 and luminosity of 1/4686 = 9 x 10 40 ergs -1 . The 
line at 3203 A is significantly fainter and its wavelength and 
linewidth are less constrained, but its presence is thought 
to confirm the identity of the 4686 A line. The ratio of line 
strengths limits the internal extinction to E(B — V) < 0.08 
mag. Surprisingly, there are no hydrogen or other emission 
lines detected in either epoch of spectroscopy. The 4686 A 
line has faded by a factor of 10 by the second epoch (at 
tpeak + 254days), while the 3203 A line is no longer visible 
then. 

The location of the event was observed at 0.2 - 10 keV 
with Chandra 315 days after the peak in the g-band light 
curve; no source (with a spectral slope of T = 2, characteris¬ 
tic of AGN) was detected above the background to an upper 
limit of Lx < 5.8 x 10 41 ergs -1 . The event was followed up 
at 5 GHz by the Very L arge Array in 2012 but not detected 
(Ivan Velzen et aDl2013h . 


1.2 Host galaxy and BH 


PSl-lOjh was coincident with the center of a host galaxy 
at 2 = 0.1696 identified in the Sloan Digital Sky Survey 
fSDSS: [Aihara e t al.l 20 1111 and U KIRT Infrared Deep Sky 
Survey rUKIDSS: lLawrence et al.ll2007l ). whose photometry 
is consistent with a stellar mass of M* , g ai = 3.6 x 10 9 M q 
dGezari et al.l[2012l ) . The spectrum of the host galaxy shows 
no narrow or broad emission lines (indicating a lack of star 
formation and AGN activity), and can be fit by a relatively 
old stellar population 1.4-5.0 Gyr old (suggesting PSl-lOjh 
was not a core-collapse supernova). 

The morphology of the galaxy is poorly constrained 
from the low-resolution images in SDSS and UKIDSS. The 
low stellar mass inferred from the colors could suggest that 
the galaxy is a disk, while the lack of star formation and 
clear 4000 A break suggest that the host may instead be 
early-type. The Magorrian relation between bulge mass and 
BH mass allows us to make some estimates (though _only 
poor) of the mass of the galaxy’s central BH tearing fc Rixl 
|200j). Assuming that the galaxy is purely a bulge of mass 
g a i = 3.6 x 10 9 Mq, the mass of the galaxy’s central BH 
would be expected to be Mbh ~ 5 x 1O 6 M0. If the galaxy 
is instead predominantly a disk whose bulge comprises only 
~ 10% of the total mass, the BH may be ~ 10 times less 
massive, Mbh ~ 5 x 10 5 Mq. At the same time, the Magor¬ 
rian relation is poorly constrained and shows high scatter 
at lo w masses, and so these estimates may not be so vali d 
('e.g. jGreene. Ho fe Barthll2008l : Ijiang, Greene fc lldl201 il l. 
In any case, the observations of the host galaxy are not in¬ 
consistent with a central BH of mass A/bh ~ 10 5 — 10 6 Mq. 


1.3 Previous studies of PSl-lOjh 


Because the spectrum of PSl-lOjh showed He II emission 
lines but no hydrogen lines, IGezari et al.l (j2012i ) proposed 
that the event was the tidal disruption of the helium- 
rich core of a red giant that had previously been stripped 
of its hydrogen-rich envelope. Various works since then 
have i nvestigate d th e event an d this interp reta tion fur¬ 
ther. iBogdanovic. Cheng fc Amaro-Seoanel (|2014f l studied 
how tidal heating of the core and orbital decay by emis¬ 
sion of gravitational waves could lead to disruption of the 
core b y a ~ 10 6 Mq BH; the ir work builds on pr evious work 
bv lMacLeod, Guillochon fe Ramirez- Rui3 (120121 ') that inves¬ 
tigated the partial disruption of r ed giant stars (see fu rther 
discussion of these papers in 451) . IGezari et al.l fl2012l ~) con¬ 
sidered in particular a red giant whose envelope had been 
stripped b y a previous interaction with a massive BH (as 
in Davies & Kind 120051 '). and they approximated its mass 
and radius from observations of a red giant in a binary sys- 
tem whose envelope had been stripped by its companion 
(jMaxted et alJ l201ll 'l . lArmiio fe de Freitas Pacheco! J2013I') 
studied PSl-lOjh using a hydrodynamic code and estimated 

the accretion rate. _ _ 

iGuillochon. Manukian fe Ramirez-Ruid (|2014|~) instead 
proposed that PSl-lOjh was the disruption of a main- 
sequence star, appealing to simulations of AGN spectra 
where He II 4686 is brighter than Ha; these authors per¬ 
forin no ph otoionizat ion calculations themselves, though 
iGaskell fe Roias Lobosl d2014ll do an d state t hat their 

findings support Guillochon. Manukian fe Ramirez-Ruid 
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d20l4) ’s idea. They claim that “over-ionization” of hydro¬ 
gen can suppress Ha emission relative to He II 4686. We 
investigate and discuss these ideas at length in Appendix 
m and explain that this idea is erroneous: instead, a high 
ionization fraction of hydrogen produces copious Balmer 
emission as ionized h ydroge n recombines in steady state. 
We thus conclude as iGezari et al. J2012l l do, that PSl- 
lOjh is likely the disruption of a helium-rich stellar core. 
iGuillochon. Manukian fe Ramirez-Ruid (120141 1 also suggest 
the presence of gas outside an accretion disk to reprocess 
and cool the radiation. 

An investigation of the accretion physics and radiation 
pro cesses that took place in PS l-lOj h has so far b een lack¬ 
ing. IGuillochon. Manukian fc Ramirez-Ruid d2014l l do per¬ 
form numerical simulations of disruption and fallback to the 
BH (though they disrupt a solar-type star rather than a 
helium core); however, their simulations lack radiation pres¬ 
sure (and accretion physics), which we argue is critical for 
understanding PSl-lOjh. 

1.4 Proposed model and plan of the paper 

In this work, we separately draw inferences from the contin¬ 
uum photometric observations, and from the spectroscopic 
data, and bring these together to propose a physical model 
for the event. 

I 11 SJ2] we use photometric observations to infer the pres¬ 
ence of an outflow and estimate its velocity, mass, and ge¬ 
ometry; then in f[3]we build on the photometric inferences to 
study the kinematics and energetics of the outflow. In we 
study the observed He II emission lines to argue that the lack 
of observed hydrogen emission essentially rules out the pos¬ 
sibility that PSl-lOjh was the disruption of a solar-type star, 
and suggest the presence of a resonance-line-driven wind. jj5] 
investigates the possibility that PSl-lOjh was the disruption 
of a helium core, and [^describes our proposed model (sum¬ 
marized in the following paragraph). We compare the prop¬ 
erties of PSl-lOjh with those of supernovae in 33 and wrap 
up with a discussion in H5] The Appendices offer detailed 
investigation of photoionization processes in PSl-lOjh. 

We argue that PSl-lOjh was the disruption of the he¬ 
lium core of a red giant star, and propose the following 
model. Only hours following disruption, stellar gas com¬ 
pleted elliptical orbits around the BH. The gas collided with 
itself and shocked, forming a highly optically thick region 
around the BH, and producing huge quantities of photons 
that were trapped by electron scattering. The majority of 
the gas and radiation fell directly into the BH, but a small 
amount of gas (~ O.OIMq) escaped at a surprisingly slow 
velocity (~ 1000 kms -1 ), likely driven by radiation pressure. 
Very little of that radiation was able to escape until the ex¬ 
panding gas became optically thin -by which time the radi¬ 
ation had already lost most of its energy to adiabatic expan¬ 
sion. We suggest that the bulk of the observed radiation in¬ 
stead had a different source: the expanding gas experienced 
another shock at larger radius that converted expansion ki¬ 
netic energy back to photons (perhaps similar to a superlu- 
min ous supernova or internal shock in a ga mma-ray burst; 
e.g.. iGal Yamll2012l : iRees fe Meszarosl[l994 1. Meanwhile, in 
the outer edge of the flow, radiation pressure trapped by res¬ 
onance lines accelerated a wisp more gas outward from the 
photosphere at a somewhat faster velocity (i.e., a line-driven 


wind). The helium recombination cascade in this wind was 
observed as broad He II lines. 

2 INFERENCES FROM THE CONTINUUM 

OBSERVATIONS, PART I 

2.1 Radius and velocity of the photosphere 

We can use the optical light curve and the temperature of 
the photosphere T e ff to estimate th e radius of the photo¬ 
sphere of the flow R p h(t) (see also IChornock et al.l 12014 : 
Ar cavi et al. 2014). We begin by considering T e g, approx¬ 
imately constrained to be 3 x 10 4 K <; T e g 5.5 x 10 4 K, as 
we explain b elow ._ 

iGezari et alJ J 2012 il use the optical and near-UV spec¬ 
tral energy distribution at f pea k —18 days and t pea k+245 days 
(where t pea k is defined as the time of peak emission in g- 
band) to infer a blackbody temperature of T e ff = 2.9 x 
10 4 K at both epochs, assuming a standard extinction law 
E(B — V) = 0.013 mag; if the extinction is greater than 
this, the temperature would be larger. These authors de¬ 
rive a second temperature constraint from the first epoch 
of spectroscopic data (t pea k — 22 days): at this time, the 
line ratio of the two observed He II lines indicates that 
the extinction is E(B — V) < 0.08 mag, implying that 
the temperature is T e ff < 5.5 x 10 4 K. IGezari et alJ ( 2012 | j 
also find that the late-time non-detection with Chandra at 
l p eak T 315 days indicates an upper limit of T e ff <; 2.5 x 10 s K 
for Lbol <j 10 44 erg s _1 . We add an additional constraint that 
T e fi ^3x 10 4 K from the presence of doubly ionized helium 
(Hc ++ ) at both epochs of spectroscopic data (f pea k — 22 days 
and t pea k + 254 days) (see IjBlIl . For simplicity, we assume 
that the temperature is constant in time, and we make cal¬ 
culations for two fiducial temperatures T e fi = 3 X 10 4 K and 
5.5 x 10 4 K. 

We now estimate the radius of the photosphere, 

Rph{t) = (^vB„(T eS (t))) ’ (1) 

where is the luminosity in a given waveband of fre¬ 
quency v, and vBv is the blackbody function. In Figure 

3 we p lot -Rp h(l) u sing the g, r, i, and z-band data from 
Gezari et al.l J2012ll . Time is plotted starting from an as¬ 
sumed zeropoint at t pea k — 67 days, which is about 4 days 
before the first detectior0. For a given temperature, calcula¬ 
tions from the different optical bands give similar values of 
R p h(t) (as expected if the dominant opacity is produced by 
a wavelength-independent process like Thomson scattering). 
A lower assumed temperature of course gives larger values 
of R p h, from i? p h(t P eak) ~6x 10 14 cm for 3 x 10 4 K, down 
to R p h(t p eak) ~ 4 x 1C) 14 cm for 5.5 x 10 4 K. 

Overplotted are curves depicting expansion at a con¬ 
stant velocity R = vot for two values of vo close to 

2 The cadence of the Medium Deep Survey is about 3 days 
llKa iscr et a l IB . and later detections take place at spacings of 
~ 3 — 12 days, so it is reasonable that the beginning of the event 
was 4 days before the first detection. If the event began slightly 
earlier, the growth of the photosphere would be slightly slower 
than linear with time; if it began slightly later, the growth would 
be slightly faster than linear with time. 
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1 10 100 1000 
time (days) 


Figure 1. Radius of the photosphere calculated from equation 
CD using the different optical photometric datapoints and as¬ 
suming a constant blackbody temperature T e ff = 3 x 10 4 K (top) 
and 5.5 x 10 4 K (bottom). Green asterisks represent g-band data, 
black diamonds r-band, orange triangles z-band, and red squares 
z-band. The black line in each panel represents expansion at a 
constant velocity, R = vq t with vq = 1100 km s —1 (top) and 
vq = 700 kms -1 (bottom); the beginning of the expansion (t = 0) 
is taken to be 67 days before peak luminosity. The blue tick marks 
denote the times when a spectrum was taken. 


1000kms 1 . We note that this expansion velocity is sim¬ 
ilar to the escape velocity at -R p h(£peak) from a massive BH, 


Vesc(r) 


3000 


( Mbh \ 1/2 
(2 x 10 5 Mq] 


6 X 10 14 cm\ 1/2 


kms 1 ; (2) 


it is similar to but somewhat lower th an t he half linewidtljf) 
of the He II 4686 A line, 4500 kms -1 dGezari et al.ll20lT 
Note that this analysis assumes that the observed emis¬ 
sion is blackbody, although the dominant opacity is due to 
Thomson scattering rather than true absorption. This effect 
may increase the size of the photosphere by a factor of a few 
above our estimates in this section. Please see Appendix [Cl 
for further discussion. 


2.2 Geometry 

The area, and hence the inferred radius, of the photosphere 
is clearly seen to rise and then fall (as a consequence of 
the rising then falling light curve, at constant temperature). 
We propose that the flow comprises an expanding narrow 
shell of dense gas enclosing a lower-density region inside. 


3 See Jb3] for discussion of the half versus full linewidth. 


The shell has mass M s h e ii, radius R B heii(t), width AR(t), 
and density p s heii(l)- We plot our proposed density profile 
schematically in the upper panel of Figure[2] Early on, while 
we see the photosphere grow, the shell is optically thick, 
and so the photosphere lies at the outer edge of the shell 
(R P h ~ Rsheii)' the ph otosp here is growing with the expand¬ 
ing shell. ( Strubbe fe Quataertl[ 2 OO 9 I referred to this situa¬ 
tion as “edge-dominated.”) This is shown by red, orange, 
and green squares in the lower panel. Later on, expansion 
causes the density in the shell to fall sufficiently that it be¬ 
comes optically thin, and so we begin to see the region inside 
it; the photosphere then moves inward with time through 
this inner region (blue and purple squares in the lower panel 
of the figure). At the peak of the light curve, the optical 
depth across the shell is unity. Since there is no sharp fea¬ 
ture in the light curve at the time of maximum luminosity, 
we infer that any drop in the density at the inner edge of 
the shell is moderate. 

We can use this model to estimate the mass and den¬ 
sity in the shell, using its opacity. In the likely densitr0 and 
temperature range for th e flow, the dominant sou rce of opac¬ 
ity is electron scattering (llglesias fc Rogerslil996l f. Assuming 
the gas is composed primarily of doubly ionized helium (see 

this opacity is K es ,He = or/2m p = 0.2 cm 2 g 1 (because 
there are two electrons for every doubly-ionized helium nu¬ 
cleus, and each ion weighs 4m p ). 

At the peak of the light curve t pea k, we have r = 
r>'es.HeLA ie ] 1 1, where S s h e ii = TThe!]/(47r/~© ie ]|) is the 


surface density (assuming spherical geometry) and I? s heii 
-R p h. The mass is then 


A/shell 


47ri? ph 

/^es,He 


■ O.O1M 0 ( 


Rph 


6 x 10 14 cm/ 


(3) 


about ~ 10% of (half) the mass of a helium core (see 
and 45.11) . 

We infer that the density in the shell is 


Pshell(f P eak) 


tCes.He AR 

8 x 10 


(4) 

(a R/R^h) Uxfo^cm ) gcm " 3 ’ 


-14 (_04_y / 
l A R/Rp h If 


corresponding to a number density n He ++ ~ 1 x 10 10 cm~ 3 . 

It will be useful later to have an estimate of the number 
density of heliun^l ions in the shell, n apec t, at the time that 
the first spectrum was taken, t spe ct = t pe a k — 22 days. In 
our picture, the density in the expanding shell decreases 
with time as n s h e ii oc oc t~ 3 (assuming that the ratio 

AR/Rshen stays fixed at ~ 0.1). This gives 


4 For temperatures T ~ 10 4 — 10 s K, electron scattering is gen¬ 
erally the dominant opacity source for number densities up to at 
least n ~ 10 12 cm -3 . At the earliest times, when the density in 
the shell was perhaps 10 3 times its value at t pe aki other sources 
of opacity may have contributed as well, but still only slightly. 

5 In : l~\1 we test ideas about the composition of the gas by as¬ 
suming that the gas has cosmic abundances, rather than being 
pure helium. In this scenario, where the bulk of the gas is hydro¬ 
gen, the mass density in the shell would be half, and the hydrogen 
number density would be twice, the values we infer assuming pure 
helium: i.e., 8 X 10 10 cm -3 or 1 X 10 11 cm -3 , for T e g = 3 x 10 4 K 
or TeS = 5.5 x 10 4 K. 
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Figure 2. Upper panel: Schematic of our proposed density pro¬ 
file for the outflow from PSl-lOjh. Curves from top to bottom rep¬ 
resent density profiles at times spaced by factors of 2. The density 
profile is composed of a lower-density power-law inside a denser 
shell of gas of fixed mass (whose width has been enhanced for 
clarity). Lower panel: Optical depth profiles corresponding to 
the density profiles in the upper panel, showing the optical depth 
from a given radius r to infinity. A horizontal line designates the 
location of the photosphere where r(> R p h) ~ 1. The colored 
squares represent the location of the photosphere as a function 
of time, from red to orange to green to blue to purple. First the 
photosphere is confined to the edge of the shell and the photo¬ 
sphere moves out with time (red, orange); then the shell becomes 
optically thin (green) and the photosphere moves in with time 
(blue, purple). 


n sp ect ~4x 10 10 cm 3 (5) 

(which we’ll take as our fiducial value) for T e ff = 3 x 10 4 K, 
or n sp ect ~6x 10 10 cm -3 for T e g = 5.5 x 10 4 K. 


3 INFERENCES FROM THE CONTINUUM 
OBSERVATIONS, PART II 

3.1 Kinematics 

We have shown that during the light curve’s rise, the growth 
of the photospheric radius is close to linear in time (with a 
velocity Vo ~ 1000 km s^ 1 ), and we hypothesized that the 
expanding photosphere is coincident with the edge of an ex¬ 
panding shell of gas. In Figure[3]we study the kinematics of 
gas particles in the gravitational potential of BHs of vary¬ 
ing mass, assuming ballistic motion. We overplot our results 


6 

o 

K 



t (days) 


Figure 3. Black curves: ballistic trajectories of particles in the 
gravitational potential of BHs of varying mass, assuming no ac¬ 
celeration due to radiation pressure. Trajectories are fixed to have 
velocity llOOkms -1 at radius 6 X 10 14 cm, to match our results 
for fiph at peak for T = 3 X 10 4 K (Figure [TJ. Dashed curve 
is for Mbh = 10 q Mq, dot-dashed curve is for Mbh = 1O 2 M0, 
and curves are logarithmically spaced by 10 1 / 2 in between. Green 
asterisks represent R p h (again from Figure [TJ. Best agreement 
is for the smallest BH masses, indicating that Mbh was small 
0$ lo3M e) or that radiation pressure was important (for weaken¬ 
ing the effective gravity of the BH on the particles). The observed 
recession of the photosphere does not at all match the shape of 
the fall exhibited by bound trajectories (Mbh )> 10 s Mg). We 
propose that the photosphere appears to recede because the gas 
density decreases, rather than corresponding to physical motion 
of particles inwards. 


for R p h from Figure |T] assuming T = 3 x 10 4 K, and fix the 
trajectories to match vo at peak radius R p h- 

The ballistic trajectories are close to the data for low 
BH mass, where the BH’s gravity has little effect on the par¬ 
ticles’ motion. This suggests for PSl-lOjh either that Mbh 
is small (<) 10 3 Mq) or that radiation pressure was impor¬ 
tant (for reducing the effective gravity of the BH on the 
particles). However, it would be surprising to find such a 
small central BH in this reasonable-sized host galaxy (e.g., 
lHaring fc Rbd 12004) ; moreover, the tidal disruption of a 
star by such a small BH would lead to incredibly super- 
Eddington feeding rates (see equation [28]), making radiation 
pressure even more important. For these reasons, we argue 
that this figure indicates that radiation pressure is impor¬ 
tant for the motion of gas particles in PSl-lOjh. 

The observed gradual recession of the photosphere after 
67 days has a very different shape from that of the plummet¬ 
ing fall exhibited by bound trajectories (Mbh )> 10 5 Mq). 
We conclude that the photosphere does not recede because 
of ballistic inward motion. Rather, we propose that after 
the shell becomes optically thin at t pea k, we are able to see 
through it outflowing gas with a run of density p(r) that 
allows us to see deeper (in the Euclidian sense) towards the 
BH as time increases. 
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3.2 Energetics 

We can construct a picture of the energy flow in PSl-lOjh 
by estimating energy scales fo r PSl-lOjh: radiatiorjf], kinetic, 
and potential. lOezari et al.l ll2012h assume a constant tem¬ 
perature of 3 x 10 4 K to estimate the energy in radiation 
that escaped from the flow, 


E ra{ ],bol, out — 


Lboi(t) dt ~ 2 x 10 51 erg. 


( 6 ) 


For comparison, we perform a similar calculation for the 
total energy emitted only in the observed bands (making no 
assumption about temperature), and find that this quantity 
is about 3 times smaller, 


E ra d,o 


Lob s {t) dt ~ 6 x 10 50 erg ; 


(7) 


since this is a firm lower limit, and in any case these quan¬ 
tities are the same order of magnitude, we take E ra d, ou t ~ 
2 x 10 51 erg to be at least approximately correct. 

We discuss how much radiation energy may have been 
present in the gas at early times in 46.2.11 The kinetic energy 
of the expanding shell of gas is 

Ek = -M a hellUo 



The kinetic energy is four orders of magnitude less than 
the radiation energy, Ek ~ 10 _4 E ra d,boi,out, very different 
from supernovae, which have Ek ~ 10 2 E ra d,boi,out dSmarttl 
l2009h . or even superluminous supernovae, which have Ek ~ 
E ra d,boi,out- We discuss this further in 46. 2. H and SJ7] We note 
that Ek oc .R 4 h , so that the inferred peak photospheric ra¬ 
dius would have to be ten times larger than our estimate for 
the energy ratio of PSl-lOjh to match that of superluminous 
supernovae. 

If we now assume that the gas was launched from 
the innermost stable circular orbit (ISCO) around a non¬ 
spinning BH (at -Ri avmc h ~ 3Rs, with Rs = 2GMbh/c 2 the 
Schwarzschild radius: see eq. m, we infer that by the time 
the gas was far from the BH, it had gained potential energy 


AEp 


CMppMsbeW 


Rlaunch 


3 x 10 51 



(9) 


similar in magnitude to the radiation energy and far ex¬ 
ceeding the kinetic energy. The total energy observed from 
PSl-lOjh is thus 


E 


tot, out 


-^rad , bol, out + Ek + AEp 

5 x 10 51 erg, 


( 10 ) 


assuming our fiducial parameters. 

Finally, supposing that PSl-lOjh was indeed the dis¬ 
ruption of a star, we estimate the energy released by the 


6 The radiation energy is much greater than the thermal energy 
of the gas, so we neglect the latter. 


accretion of (half of) a stellar mass of gas onto a BH, 


Edis 


r/MtC 2 


2 



( 11 ) 

( 12 ) 


where the factor of 1/2 accounts for only the bound half of 
the disrupted star, rj ~ 0.1 is a typical efficiency of convert¬ 
ing accretion into radiation, and we have normalized to the 
mass of the helium core of a solar-mass red giant star (see 

SO- 

Thus we uncover a hierarchy of energy scales, 


Ek <C E ra d ,bol,out ~ AEp < Edi s . 


(13) 


These energy scales offer insights and puzzles for our under¬ 
standing of PSl-lOjh, which we describe in 1 )6.21 


INFERENCES FROM THE HE II 
OBSERVATIONS 


4.1 Not a main-sequence stellar disruption 

As described in the introduction, the optical spectrum 
of PSl-lOjh shows the emission lines He II 4686 A and 
3203 A, but_ no hydrog en lines—in particular, no Ha or 
H/3. iGezari et al.l (120121 1 interpreted these observations to 
mean that PSl-lOjh was the disruption of a helium- 
rich ob ject, su ch as the stripped core of a red gi¬ 
ant star. Guillochon, Manukian fc Ramirez-Ruij (l2014ll and 


iGaskell fe Roias Lobosl ( 2014i ) instead propose that PSl- 


lOjh was the disruption of a main-sequence star, based on 
photoionization calcul ations for AGN (u sing the publicly 
available code Cloudy: iFerland et al.ll2013t ) ; they argue that 
the observed line ratio lower limit of He II 4686 to Balmer 
Ha can occur for photoionized gas of cosmic abundances, 
and does not require an enhanced abundance of helium. 

In Appendix El we extensively investigate these claims 
by performing and analyzing our own Cloudy photoion¬ 
ization calculations. Here we briefly outline the important 
physical processes and then discuss results. The gas where 
the lines are produced has high enough density to be eas¬ 
ily in photoionization equilibrium: every photoionization is 
quickly balanced by a recombination, in which the ion joins 
an electron, and the electron “cascades” down to the ground 
state, releasing line-transition photons as it goes. This mech¬ 
anism by far dominates line emission of He II 4686 and Ha; 
the idea that the hydrogen could be “over-ionized” and thus 
produce less emission is incorrect—more ionized gas pro¬ 
duces more recombination events and so more line emission. 

One way to suppress Ha relative to He II 4686 is to have 
a signihcant population of neutral hydrogen that absorbs Ha 
photons—but few enough (singly-ionized) helium ions that 
He II 4686 is not too absorbed as well. Relatively high gas 
densities and large incident fluxes generally give rise to the 
most relevant conditions. 

However, our Cloudy calculations uncover significant is¬ 
sues, especially: 

(i) Th e cal culat ions referenced a nd prese nted by 
jGffiUochon J _Manukian_& Ramirez-Ruizi (|2014| 'I and 
iGaskell fe Roias Lobosl (2014) are for stationary gas, 
but the broad linewidth of He II 4686 indicates large 
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velocity gradients, which would diminish optical depth 
significantly: Ha would not be effectively absorbed and 
suppressed; 

(ii) Even if we make the same assumptions as these au¬ 
thors, we find that the ratio of He II 4686 to Ha doesn’t 
reach the observed lower limit for any density and incident 
flux level, especially including the appropriate parameters 
inferred from the continuum observations m-, 

(iii) Most problematically, Cloudy has serious conver¬ 
gence issues in much of the region of parameter space where 
these authors claim PSl-lOjh would lie (high density and 
large incident flux), and so the numerical results likely are 
not trustworthy anyway (R. Porter, personal communica¬ 
tion). We discuss this further in Appendix [A] 

In view of these results, we argue that it is highly un¬ 
likely that the line-emitting gas of PSl-lOjh had the cosmic 
abundances of hydrogen and helium: thus it was probably 
not the disruption of a solar-type star. Instead we argue for 
the disruption of a helium-rich obj ect such as t h e strip ped 
core of a red giant, as suggested bv lGezari et al.l (120121 '). 


4.2 Origin of the He II 4686 A line 

The He II 4686A line is produced through the recombination 
cascade (n = 4 to n = 3) when He ++ recombines with an 
electron to produce He + . Provided a region of helium gas is 
(1) predominantly doubly ionized, (2) has had time to reach 
ionization equilibrium, and (3) is optically thin to the 4686 A 
photons (see j|B]for justification of these three assumptions), 
the line luminosity will be 

T 4686 ~ {hU4eS6)oi4,6S6(EM) , (14) 

where hu^ese = 4.24 x 1CP 12 erg is the energy of the tran¬ 
sition, 04686 = 8.20 x 10 _14 cm 3 s^ 1 is the eff ective Case 
B recombination coefficienlQ at T = 4 x 10 4 K (iOsterbrockl 
1 19891 . p. 81), and {EM) = n e n He ++ V is the emission mea¬ 
sure of the line, with n e and n He ++ the number densities of 
free electrons and He ++ , and V the volume of the emitting 
region. Using the luminosity L 4686 = 9 x 10 4O ergs _1 ob¬ 
served at fpeak — 22 days, we infer that the emission measure 
of the gas producing the He II 4686 A line in PSl-lOjh is 

(EM) - 3 x 10 65 cm -3 . (15) 

We now compare this observed {EM) with the emis¬ 
sion measure of the region just outside the Thomson photo¬ 
sphere, R p h- We make the assumption that the radial thick¬ 
ness of the line-emitting region is similar to R p h, because 
of the observation that the He II lines are centered on the 
galaxy’s lines (if the emitting region instead were thin, most 
of the receding part of the flow would be blocked by the 
Thomson photosphere); the volume of the region is thus 

V ~ (few x Rph) 3 ~ 10 45 cm 3 , (16) 

' We make calculations using this temperature which is close 
to the temperature inferred for the photosphere T > 3 x 10 4 K. 
Although the gas outside the Thomson photosphere at R p h may 
not have exactly this temperature (since the gas is optically thin 
in the continuum), its temperature likely is not too different from 
this since the gas remains coupled to the radiation through line 
transitions. The recombi nation coefficient scales roughly as T~°' S 
llVerner fe Ferlandlll996h . 


where we have used R p h (Aspect) ~ 4 x 10 14 cm at t sp ect = 
Ipeak — 22 days. If the density in the gas is similar to the 
density we infer for the shell at that time (n He ++ ~ 4 x 
10 10 cm -3 ; see eq.0), the emission measure in that gas would 
be 

{EM) ~ 10 66 cm -3 , (17) 

quite similar to the emission measure inferred from the line 
luminosity. This lends additional support to our interpre¬ 
tation that the line-emitting region lies just outside the 
electron-scattering photosphere. 


4.3 Line-driven wind and ultraviolet absorption 
lines 

Broad-absorption-line quasars (BAL QSOs) and O-stars 
both show evidence for line-driven winds, in which gas at 
the outskirts of these objects absorbs luminous ultraviolet 
radiation via resonance line transitions (from ground state 
to higher energy level), and is accelerated outward by ra¬ 
diation pressure. Observational evidence for the presence 
of these winds is provided by spectra showing blueshifted 
resonant absorption lines. The spectral energy distribu¬ 
tions of O-stars are essentially those of blackbodies with 
T e fF ~ 20, 000 — 50, 000K; those of BAL QSOs show a simi¬ 
lar blackbody-like UV bump, supplemented by strong X-ray 
emission. By analogy to these objects, we expect UV-bright 
PSl-lOjh to have produced a line-driven wind as well. We 
make some estimates about this wind and offer some pre¬ 
dictions for observations of similar events in the future. 

In H2.ll we inferred the presence of a shell, whose ex¬ 
pansion (at Vsheii ~ 1000 km s -1 ) we suggest is driven by 
radiation pressure from photons trapped by Thomson scat¬ 
tering. In the outer reaches of this shell, beyond R p h, the 
gas is optically thin to Thomson scattering. This outer gas 
should be subject to the same line-driven instability that 
drives high-velocity winds from O-stars. Since the observed 
spectrum shows us that He + was present in PSl-lOjh, singly 
ionized helium is likely to provide a substantial line force 
for driving the wind; in HB2I we show that the He Lyman a 
transition at Eueha = 40.8 eV is highly optically thick. By 
analogy with O-stars, we expect that there are many reso¬ 
nance metal lines that contribute to driving the wind; the 
resultant accelerations ensure a large radial velocity gradi¬ 
ent in the gas, broadening the lines substantially. In ' ll S3 1 we 
speculate about the mass-loss rate of the wind. 

Line-driven winds from O-stars have terminal 
velocities somewhat larger than the escape veloc¬ 
ity at the win d’s launc h radius, typic ally by a fac¬ 
tor of 3 llCastor. Abbott fe Klein 19751: Abbottl Il982l : 


iHowarth fe Prinial Il989l : iLamers fe Cassinellil Il999f) : radia¬ 
tion pressure gives the gas order unity of its gravitational 
binding energy Th e same is be lie ved to be true for BAL 
QSO winds (IMurrav et al.l Il995l : iProea. Stone fe Kallmanl 
120001) . In PSl-lOjh, the observed half-width of the He II 
4686 line corresponds to viinewind ~ 4500 km s~ 4 , which we 
infer to be the terminal velocity of a wind. This is a factor 
of a few times faster than the inferred expansion velocity 
of the continuum-driven shell. We thus suggest a picture 
in which the outer layers of gas in the shell (just beyond 
Rph) are additionally accelerated by radiation pressure on 
ultraviolet resonance lines to velocities somewhat faster 
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than the shell’s velocity. The observed He II emission 
lines are likely produced in the outer parts of the shell 
(as suggested in S2]) and in the faster (lower-density) 
line-driven wind beyond. The observed velocity of the 
line-emitting gas together with the size of the launch radius 
(~ R p h) offers another estimate of the mass of the BH, 

, r ^linewind^ph 

m bh ~ - — - 

~ 3x 11)51,0 (^rrfeir) (18) 

This inferred BH mass is remarkably similar to the BH mass 
we infer in asm via the independent argument of requiring 
that the disruption of a helium core take place outside the 
BH’s innermost stable circular orbit. 

We note that the half line-width (rather than the full 
line-width) is an appropriate measure of the wind velocity 
for the following reason. In an optically thin outflow, we 
observe emission from gas moving towards us (+uii ne wind) 
and away from us (—vunewind), which together produce a 
broad emission line centered on the velocity of the central 
BH, whose width corresponds to twice the outflow velocity. 
We find that this is the situation for the He II 4686 line of 
PSl-lOjh, be cause the line is centered on the galaxy’s lines 
llGezari et al.ll20l1 l. and because the gas is marginally opti¬ 
cally thin to He II 4686 photons (see Appendix 0. If there 
were significant 4686 A absorption in addition to emission, 
the situation would be different: any observable absorption 
would take place along our line of sight, i.e., only in gas 
moving towards us; this would produce a P-Cygni line pro¬ 
file, which comprises blueshifted and redshifted emission su¬ 
perposed with blueshifted absorption (generally appearing 
simply as blueshifted absorption and redshifted emission). 
However, there was not enough He + gas in the n = 3 en¬ 
ergy level to produce significant 4686 A absorption because 
the gas was too cool (see Appendix 0; thus the He II 4686 
appeared as a pure emission line. 

We expect that if an early-time ultraviolet spectrum of 
PSl-lOjh had been taken, the He La line would have been 
observed. In contrast to the He II 4686 line, the He La line 
would have exhibited a P-Cygni profile, because there was 
much more Hc + gas in the ground state n = 1 than at n = 3, 
and so the gas was highly optically thick to He La absorption 
(see Appendix0. Absorption of He La takes place between 
us and R p h, through gas moving towards us with veloci¬ 
ties ranging from ~ 0 — uiinewind! hence the absorption line 
would appear broad and blueshifted by ~ 4500kms _1 . The 
gas outside _R p h all emits He La as well, so a broad emission 
line (centered on the galaxy’s lines) maybe have been visible 
superposed with the blueshifted absorption. More specula¬ 
tively, perhaps the width of the line (prior to f pea k) decreased 
with tim(0 as An ~ v eac (R v h) oc R ~ h 1//2 oc t -1 / 2 . 

Depending on the composition and ionization state of 


8 Note that the observations of the He II 4686 li ne at ^ + 
254 days were unable to measure its late-time width llGezari et alj 
1 HH).’ We speculate that the width of this line (and any other 
putative lines, such as He La:), may have followed the velocity of 
the photosphere (perhaps decreasing prior to t pea k and increasing 
afterwards), but future detailed study of velocity evolution in line- 
driven TDE outflows will require more sophisticated models. 


the gas, other ultraviolet resonance lines would have been 
visible as well, also with P-Cygni profiles (again, because 
they are due to absorption f rom the ground state n — 1 ): 
e.g., C III, Si IV, C IV. fSee IStrubbe fc QuataertJboill for 
detailed spectroscopic predictions of tidal disruption events 
in a slightly different context.) Other Lyman series transi¬ 
tions of He ++ may also have been visible, and perhaps even 
a weak hydrogen Lya line, if a trace amount of hydrogen 
was present. PSl-llaf, a more recently discovered tidal dis¬ 
ruption candidate, was observed with near-UV spectroscopy, 
and exhibited two broad and deep (unidentified) absorption 
lines at ~ 2000 — 3000 A 24 d ays after the event was dis¬ 
covered dChornock et al.l 120141 4 ; the physics was probably 
similar in this event. Efforts are underway to obtain ultra- 
violet spectroscopy of disruption candidate ASASSN-14ae 
dHoIoien et al.l2014h using the Hubble Space Telescope (S.B. 
Cenko et ah, personal communication). 


5 TIDAL DISRUPTION OF A HELIUM CORE 

We have laid out evidence that PSl-lOjh was too helium-rich 
to be the tidal disruption of a solar-type star. As we explain 
further in 0 the high fraction of energy output in radiation 
relative to kinetic energy (i m, along with the behavior of 
the temperature and radius of the photosphere f fj 2 .ll) . are 
unlike the behavior of any known type of supernova. Now 
we investigate the interpretation that PSl-lOjh is the tidal 
disruption of the heliu m c ore of a red giant star (proposed 
first bv lGezari et al.ll2012l h 

5.1 Tidal disruption radius and black hole mass 

We calculated red giant models using the open-source stel¬ 
lar evolution code Modules for E xperiments in Stellar As¬ 
trophysics (version 4219 of MESA. [Paxton et al.ll201lh . Our 
fiducial model has initial mass IMq, initial metallicity z = 
0.02, and age f age = 1.24 x 10 10 yr. Red giant cores are 
very dense: the fiducial model has radius R co re ~ 0.03 -Rq 
and mass M cole ~ 0.3AI 0 , giving a mean density p c ore ~ 
2 x 10 4 gem -3 . 

Because the core’s density is so high, the tidal force 
required to disrupt it is much greater than that for a solar- 
type star. We estimate the core’s tidal disruption radius (i.e., 
the distance from the BH at which the BH’s tidal gravity 
can fully disrupt the corefj 



O If ) 11 ( ^ core ) ( Mcore A ^ f Mbh \ ^ 
X (o.O3R 0 ) \0.3M e ) (2x10 5 Mq) 

9 The tidal disruption radius of a star is R t = ^(Mbh/M*) 1 ^ 3 
where £ is typi cally order unity. Simulations by 
iGuillochon Sz Ramirez-Ruid d2013h of disrupting solar-type 
stars indicate that (depending somewhat on stellar structure) 
full disruption takes place for £ « 1; we make the assumption 
that this result holds for red giant c ores a s well. Note that 
I MacLeod, Guillochon &; Ramirez- Ruid (|2012|V s work on the 
disruption of red giant stars studies the disruption of red 
giants’ low-density envelopes rather than the dense cores we are 
interested in here. 


cm 


© 2015 RAS, MNRAS 000. mf26l 

































Tidal Disruption Candidate PSl-lOjh 9 


/ Rears \ f M core \ 1/3 / Mbh \ 

\0.03R Q ) \0.3M Q J ^2x10 5 Mq) 


where Rs = 2GMbh /c 2 is the BH’s Schwarzschild radius. 

The requirement that the core be fully disrupted to 
produce an observable event places significant constraints 
on the mass of the BH. If the BH were too massiv<F°l 
Mbh ;> 1 O 6 M 0 , the core’s tidal radius would lie inside the 
BH’s event horizon (_Rt, core ^ Rs) and the BH would swal¬ 
low the core whole rather than tidally disrupting it. If the 
star’s orbital pericenter lies outside the event horizon but in¬ 
side the BH’s ISCO (risco = 3Rs for a non-spinning BH), 
much of the gas may be acc reted by the BH shortly af¬ 
ter disruption (see [Haas et a] .1120 1 2 for simulations of white 
dwarf disruption at different pericenter distances). It is likely 
that this scenario would not produce much emission either, 
although future work incorporating general relativistic ef¬ 
fects, accretion physics, and radiation will be necessary to 
make firmer predictions. Based on these arguments, we infer 
that Mbh ^ 10 6 Mq —and more likely Mbh ~ 2 x 10 5 Mq 
or less. We adopt this value of Mbh = 2 x 1O 5 M 0 as our 
fiducial BH mass, with a fiducial tidal disruption radius of 
-ftT.core ~ 2 X 10“ Cm. 

The inferred bolometric luminosity Lbol ^ 2 x 

10 44 erg s ” 1 ( iGezari et al.l 120121 1 is thus about 4 times the 
Eddington luminosity appropriate for doubly-ionized helium 
gas, 


TEdd,He 


4ttGMbhc 


^es,He 

5 x 10 43 


f Mbh \ 
^2 x 10 5 A4© J 


ergs 1 . 


( 20 ) 


This indicates that radiation pressure plays a significant role 
in the dynamics of PSl-lOjh, as we shall see further shortly. 


5.2 Fallback 


In the_ now _((classical” p icture (e.g., 

lLacv. Townes fc Hollenbacbl 1 19821 ; iReed 1 1988T ) . as the 
stellar core disrupts, half of its gas becomes bound 
to the BH. The bound gas departs the BH on highly 
eccentric elliptical orbits. The most bound gaE] has 
energy £_ ~ (3GMb h / Rt. cor e )(R core /Rt, core) 

Lac y. Townes fc Hollenbacbl 1 19821 ; iLi, Naravan fe Menoul 
200211 . so that its apocenter is 


Ra 


RA 


3R C 


6 x 10 1 


( 21 ) 


Rt,c 


,2 x 10 11 cm 


\ 2 / Rcore 

) \0.03R q 


(roughly twice the semi-major axis because the orbit is 
highly eccentric). This most tightly bound gas begins re- 


10 This upper mass limit would presu mabl y be som ewhat 
larger i f we accou nt for spinning BHs; s ee iKesdenl 120121 . Work 
by iBogdanovic, Cheng &; Amaro-Seoanel (l2014h also suggests a 
weaker upper limit to Mbh by allowing tidal heating of the core 
to reduce the core’s density; see 

11 We assume for simplicity that the star is maximally spun up 
during disruption. 


2 ^ 3 turning to pericenter after the “fallback time” 

,(19) 

' I?T,core \ 3/,a f Rt,c ^ 


^fallback 


2 tt 

63G 


/ Rt.cots \ 
V R core / 


V GMbh ) 


7 x 10 J 


( Mbh ^ 

1 ( -^T,core ^ 

l 2 x 10 s Mq J 

1 V 2 x 10 11 cm/ 


Ra 


- 3/2 


O.O3f ? 0 1 S ’ ^ 

(e.g., iR.eesi fliissl : ILi. Naravan fc Menoul [2002 1—only about 
2 hours. Bound gas continues “falling back” to the BH at a 
rate 

1 M c ,— / t \- 5 / 3 


Alfa,! [ back ~ “ 


( — ) 

V ^fallback / 


(23) 


3 ^fallback V ^fallback / 

<jRee jjl 9881: IPhinnevI fi989l R depicted by the solid black 
curvJ 43 ! in Figure |4] Note that most of the mass falls back 
at early times: within a time t ~ 10tf a iiback ~ 1 day after 
disruption, 80% of the bound gas has fallen back to pericen¬ 
ter; already half of it has fallen back after only t ~ 3£f a iiback- 
We argue in 1 )6. II that this could be an explanation for the 
shell-like geometry of the expanding gas. 

When returning gas reaches the vicinity_of the 

BH, it is expecte d to sh ock on itself 


f (e.g., _ 

Evans & Kochanck 1989; Kochanck 1994; Guillochon et al. 


Rees 


2009ll and convert orbital kinetic energy into radiation en¬ 


ergy with some efficiency r/. If this radiation were able to 
escape immediately, the bolometric light curve of the event 
would be proportional to the mass fallback rate Mf a nback- 
The peak luminosity of the event would be 

2 


A ( /. f;i,| I back ) 


r/ Lfal I back ' 


10 4 


Mbh 


0.1 A/fallback (^fallback ) C 
1/2 


(24) 


2 x 10 s M 0 

f R'l/core A 

V 2 x 10 11 cm/ 


with an effective temperature of 

0■ 1 A/fallback (^fallback ) C 


M c , 


O.3M 0 


R a 


0.03 R P 


3/2 


ergs 


kT e ft 


k 


^^l.core^SB 


Mbh 


2 x 10 5 Mq 


\ 1/8 / M cole \ 1/4 / 

• J \ O.3M 0 J { 


Rt,, 


2 x 10 11 cm/ 


Re, 


O.O3R 0 


3/8 


keV 


(assuming the radiation had time to thermalize), lasting for 
a timescale of ~ few x ffaiiback ~ hours. Here ctsb is the 
Stefan-Boltzmann constant, and Ffaiiback is the rate at which 


12 There may be slight or moderate deviations from this canon- 
i cal form depending on details of t he core’s d ensity profile: see 
iLodato. King fc Pringl3 ll200flh and [Guillochon & Ramirez-Ruiz 
(1201311 for results for solar-type stars. 

13 In Figure |4] we plot the fallback rate from equation (f23l) 
multiplied by a factor 3/r(l/3) x exp[— (tfaiiback/^) 2 ] to ^p- 
proximate the onset of f allback; this shape is similar that of 
I Evans &; Kochanekl <|l989h ’s numerical result for the disruption 
of a solar-type star. This diminishes the maximum fallback rate 
by a factor of several. 


- 5/4 

(25) 
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t (days) 


Figure 4. Black solid curve: Rate of energy generation in gas 
falling back to pericenter in units of the Eddington luminos¬ 
ity, 0.1Mf a nback c2 /^Edd? for our fiducial model of a red gi¬ 
ant core (Mcore = 0.3 Mq and R c ore = O.O 3 R 0 ), assuming 
Mbh = 2 x 10 5 M o and RT,core = 2 x 10 11 cm. The fallback 
rate is hugely super-Eddington shortly after disruption, and re¬ 
mains super- Eddington for mont hs. Green diamonds: g-band data 
for PSl-lOjh dGezari et al.11201211 in units of the Eddington lumi¬ 
nosity. That the light curve peaks ~ 3 orders of magnitude later 
than the fallback rate peaks indicates that the photons generated 
in the fallback shock get trapped (though more photons may be 
generated in a shock at larger radius; see £16.2.11) . Photons can¬ 
not escape until much later when the flow becomes optically thin 
(at t ~ 67days). After that time, the light curve approximately 
follows the energy generation rate. 


orbital kinetic energy is dissipated in the shock at pericen- 
ter at time tfaiiwk, assumed to be ~ 0.1Mf a iib ac k(if a iib a ck)c 2 
for some reasonable value of the efficiency factor > 7 . This 
timescale, luminosity, and temperature are dramatically dif¬ 
ferent from those observed for PSl-l Ojh. _ 

Recent numeric al work (e.g., IShiokawa et al.l 120151 : 
iBonnerot et al.ll2015l l suggests modifications to this “clas¬ 
sical” picture: the shock at pericenter may be rather weak 
compared to stream-crossing shocks that occur near apoc- 
enter (i? apo : see eq. 1211 . The effect is to prolong the dura¬ 
tion of peak accretion, by a factor of 3-10, and to reduce 
the peak accretion rate. However, even with the prolonged 
timescale, the accretion timescale is still two orders of mag¬ 
nitude shorter than the observed duration of peak lumi¬ 
nosity. Additionally, the shock at pericenter may be strong 
nevertheless for disruptions close to the ISCO, as PSl-lOjh 
seems likely to be. We continue here assuming fallback times 
and accretion rates that correspond to a shock at pericenter 
rather than apocenter, but discuss these issues further in 
38^21 

5.3 Trapped radiation 

In Figured we plot the g-band light curve observed for PSl- 
lOjh (green diamonds) over the classical estimate for the 
energy generation rate at the shock, ~ 0.1Mf a iib a ckC 2 (see 
Footnote m- The different shapes are striking: the light 


curve peaks ~ 3 orders of magnitude later in time than 
the energy generation rate, and the bolometric light curve’s 
peak brightness is ~ 3 orders of magnitude below the peak 
energy generation rate (assuming that the g-band luminos¬ 
ity is ~ 10% of the bolometric luminosity). Interestingly, 
the light curve’s decay does approximately follow the en¬ 
ergy generation rate. 

We thus propose that, instead of immediately escaping, 
most of the radiation produced either in the shock at peri¬ 
center or by the accretion disk was trapped at small radii by 
electron scattering. The optical depth in the flow at pericen¬ 
ter would initially have been enormous, 


r(3t ia 


^es,He (-Mcore/4) 

47 r ^T,core 


6 X 10' 


f Mcore j 

{ -f^TjCore \ 

V 0.3M© J 

V 2 x 10 11 cm/ 


(26) 


> 1 , 


(where we are approximating the flow as spherical). If the 
stellar gas around the BH remained static, the time for pho¬ 
tons to diffuse out would be 


tdiff(RT ,core) 


T ( 31 fallback ) -^T ,c 


(27) 


10 


/ Mcore \ 

f -^T,core \ 

l 0.3M© ) 

V 2 x 10 11 cm/ 


yr, 


much longer than the fallback time or dynamical time (or 
rise time of the observed light curve). 

We propose that the trapped radiation pushes gas out¬ 
ward: the radiation produced via fallback is easily energetic 
enough to unbind a fraction of the originally bound half of 
the stellar gas, since the fallback rate shortly after fallback 
begins (say, at time 3tf a iib ac k) is orders of magnitude above 
the Eddington rate (see Figure [4]), 

Mfallback ^ g ^q 3 f Mbh \ ^ / M co re \ 

Me dd \2 x 10 s Mq J (o.3Mq J 

/ RT.core \ ~ 3 / _gcor^\ ' ( } 

V 2 x 10 11 cm/ \0.03R q J v ' 

The observations directly show that there is enough energy 
available to lift a substantial amount of mass out of the 
potential well of the black hole; the energy we observe in ra¬ 
diation is sufficient to lift ~ 0.01M© from a radius of several 
Rs (c.f. eqs. |5]and [SJ). 


6 INTERPRETATION 

We now draw together the arguments from the previous sec¬ 
tions to propose the following picture for what took place in 
PSl-lOjh. 

6.1 Simple model 

The helium core of a red giant star was disrupted by a rel¬ 
atively low-mass BH (Mbh $ 2 x 10 s Mq). In the hours 
following disruption, newly-bound stellar gas completed el¬ 
liptical orbits around the BH and fell back to pericenter at a 
highly super-Eddington rate (eq. 1281 , shocking on itself and 
quickly generating huge amounts of radiation. The gas was 
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incredibly optically thick lea. 1261) . and so essentially all of 
the radiation was trapped in the gas by electron scattering. 

The radiation contained enough energy to gravitation¬ 
ally unbind a significant fraction of the stellar gas that had 
fallen back to the BH: about 10% (O.OIMq ~ 0.1(M core /2)) 
was driven out in an optically thick shell at v ~ 1000 kms^ 1 , 
whose expanding edge was observed as the expanding pho¬ 
tosphere -Rph- While the shell remained optically thick, pho¬ 
tons were advected out with the gas faster than they could 
make progress diffusing out. Radiation was mostly only able 
to escape (and reach our telescopes) once the shell had ex¬ 
panded enough that the density had fallen enough to allow 
the shell to transition from optically thick to optically thin. 
This transition took place around t ~ 67 days after disrup¬ 
tion, producing the observed peak of the light curve (long 
after the peak of accretion). 

Outside the photosphere for Thomson scattering, a frac¬ 
tion of the radiation was trapped by resonance lines (e.g., 
helium Lyman a), which accelerated a tiny wisp more gas 
to velocities ~ 4500 km s -1 . This outermost gas was pho- 
toionized by the radiation coming from the Thomson pho¬ 
tosphere, and produced the broad He II recombination lines 
that were observed (which serve as clues to the composition 
of the disrupted object). 

Eventually the expanding shell became optically thin, 
and from then on we could see through to gas inside it. 
Although gas particles continued to expand outwards, the 
location of the photosphere now moved inwards as the den¬ 
sity of the inner gas fell and allowed us to see deeper and 
deeper inside: see m\ and Figure 0 We hypothesize that 
the shell geometry was present because of the very short 
timescale (At ~ few x finback <S t pea k) over which most of 
the mass fell back to pericenter and most of the radiation 
was generated. The gas interior to the shell fell back to peri¬ 
center at later times t few x tfaiiback: it was subsequently 
driven out by radiation pressure, and we speculate that it 
may have had more of a wind geometry p oc r -2 , as proposed 
bv lStrubbe fc Quataertl (12009 '). 

Now that the total column from the BH out through the 
flow was so much smaller, most of the radiation generated 
close to the BH could escape, rather than being swallowed. 
The radiation continued to be reprocessed by the outer gas 
(hence its cool 3 x 10 4 K temperature), but not nearly as 
heavily as while the shell was optically thick. Thus we see 
that the photometric measurements at these late times are 
similar to 0.1M fa iib ac kC 2 (Figure HR 


6.2 Energy comparisons 

This simple model explains some of the energy scales we 
inferred in T3.2I 

(i) As the shell of stellar gas expanded, trapped radia¬ 
tion energy was converted to kinetic and potential energy 
of the gas particles. Conserving energy, we’d expect to find 


14 We should more properly compare the bolometric luminosity 
Lboi with 0.1Mf a iib ac kC 2 , but as the temperature remains roughly 
constant, the 3 -band emission remains a constant fraction ~ 0.1 

°f Tbol- 


that the initial radiation content of the shell (before it ex¬ 
panded) -E ra d,o is equal to the total energy in the gas when 
the remaining radiation was finally able to escape, 

E ra d,0 r '' J E rad ,bol,out + Ek + AEp , (29) 

Our observation that E ra d,boi,out ~ Ek + AEp suggests 
that the initial radiation content of the shell was a few 
times the quantity of radiation that escaped (E ra d,o ~ 
few x E rad ,boi,out). We discuss this further in 1 16.2.11 

(ii) The total energy released by the disruption Edi s was 
almost an order of magnitude greater than the total observed 
energy of the shell (eq. flOl) , 

Etot.out = (E ra d,boi,out + Ek + AEp) ~ 0.2Edi S , (30) 

assuming that 77 ~ 10 % of the accretion energy was con¬ 
verted to radiation (see eq. HU): Evidently, the remaining 
~ 80% of the radiation energy liberated by the disruption 
was quickly swallowed by the BH, presumably entrained in 
the ~ 90% of the gas that was not blown away in the shell. 
(Note that we infer that the fraction of gas that was blown 
away in the shell is ~ 0.01ALq/0.15Mq ~ 7%, discussed 
further in H51 ) 

Simulations of accretion physics often measure a quantity 
we’ll call rjesc , the efficiency with which mass fed to the BH 
is converted to radiation that then manages to escape from 
the flow to infinity. This quantity differs from what we’ve 
called r/, the efficiency with which mass fed to the BH is 
converted to any radiation—much of that radiation can then 
be swallowed by the BH and not contribute to rfesc- Our 
estimate for r ] e sc in PSl-lOjh, from comparing the observed 
quantity Etot.out to (half) the mass of the stellar core (i.e., 
Edis /77 in eq. [TT|), is 

'■ (31) 

This result is intriguingly similar to results of the recent 
3-di mensional radiation magne tohydrodynamic simulations 
by I Jiang, Stone fe Davisl (120141) who find r/ e ac = 0.045, dis¬ 
cussed further in U5] 

(iii) In supernovae, the majority of the radiation energy is 
typically converted to kinetic energy, le ading to the ob served 
E ra d,out ~ 10 49 erg « E K ~ 10 51 erg (ISmarttl 12009 ). PSl- 
lOjh contrasts starkly with this: Ek ~ 10 47 erg -C E ra d, ou t ~ 
10 51 erg ( 93.21) . We suggest that the kinetic energy was much 
smaller than the radiation energy for two reasons. First, in 
contrast to supernovae, much of the energy initially stored 
in the radiation was converted to potential energy instead of 
kinetic energy—photons had to do work against the gravi¬ 
tational potential of a massive BH, work the photons don’t 
have to do in a supernova. Another reason for Ek -C E ra d,out 
could be that the shell experienced a shock at large radius, 
converting a substantial fraction of the kinetic energy into 
photons; see 96.2.1 1 This difference in energetics provides 
strong evidence that PSl-lOjh is not a type of supernova. 
We offer a fuller comparison with supernovae later in 

(iv) It seems significant that the outflowing shell of gas 
was only barely unbound from the BH, as we deduce in two 
ways: (1) The kinetic energy is tiny compared with the po¬ 
tential energy (Ek -C AEp: 113.21) . and (2) The outflow veloc¬ 
ity is similar to the escape velocity at large distance from our 
fiducial BH (eq.EJ). By contrast, winds from AGN or O-stars 
typically have velocities several times the escape velocity 
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from the radius where they are la unched ('e.g. jMurrav et all 
Il995l : lLamers fe Cassineli3 1 199911 : if PSl-lOjh had followed 
suit, its expansion velocity would have been ~ 30 times 
greater. Did some process regulate radiation pressure to can¬ 
cel gravity almost exactly? 

(v) The model does not explain the observed constant T e ff 
of the photosphere. The large size of the photosphere relative 
to the tidal radius f?T.core suggests that the radiation cools 
by scattering and/or advection through optically thick gas 
between the location where the radiation is generated (the 
shock or inner regions of the accretion disk) and the location 
where the radiation begins to free-stream. But as we point 
out in equation (1321) , adiabatic expansion would suggest that 
the photosphere cooled as the gas expanded, which was not 
observed. In H6.2.1I and fJS] we discuss further ideas about 
what sets the temperature. 


6.2.1 Adiabatic cooling? - An additional shock 


The final energy comparison of this section addresses the 
question of adiabatic cooling of the radiation. Most of the 
radiation escaped close to the time t pea k, when the shell be¬ 
came optically thin, so E r ad ,out is approximately the radia¬ 
tion content of the flow at that time: -E ra d(t P e a k) ~ E rad ,out- 
We have argued that at times earlier than this, most of 
the radiation was trapped by Thomson scattering within 
the expanding flow. As the flow expands, we expect the 
energ y in radiation to decrease due to adiabatic coolin g 
(e.g., Istrubbe fe Quataertl 120091 ; iRossi fe Begelmanl l2009h . 
The temperature of an adiabatically expanding photon gas 
scales as T oc R -1 ^ 3 oc R~^ ge , where V and R e dge are the vol¬ 
ume and radius of the flow; the energy in radiation therefore 
scales as 


-Erad ~ aT*V OC R ed g e 


(32) 


We can use this relation to work our way backwards in time 
and infer the radiation content of the gas early on. Supposing 
the flow originated close to the tidal radius for a disrupting 
red giant core ~ 2fTr,core fse e 115.11) . we’d estimate that the 
radiation energy was initialhl 15 l a whopping 


-Er a d,0 ~ 3x10 


f -^ph(^peak) j 

( -^T,core ^ 

I 6 x 10 14 cm ) 

V 2 x 10 11 cm) 


erg , (33) 


fully 100 times larger than the energy released in the dis¬ 
ruption of a red giant core (eq. 1111) . Where could all that 
energy have come from? 

Disrupting a much more massive star could solve the 
problem, as i?dis would be correspondingly larger; however, 
postulating the disruption of an exceedingly rare ~ 30 Mq 
core seems far-fetched. Instead, we argue that the expansion 
was not adiabatic. 

For comparison, in a typical supernova explosion, ex¬ 
pansion is adiabatic and homologous: concentric shells ex¬ 
pand at different velocities such that they do not interact 
with each other, and most of the radiation content of the 
gas is converted to kinetic energy (so that _E ra a -C -Ek ) 
as the flow expands in an orderly fashion dArn e ttl ll982h . 


15 Even using equation to estimate the radiation energy 
at the time of first detection t\ gives an uncomfortably high 
E ra d(fl) ~ 2 x 10 52 erg. 


The class of superluminous supernovae, however, which have 
_ E ra d ~ 10 51 er g, app ear to require non-adiabatic behavior 
dQuimbv et al.lfioill l. 

In contrast to the standard supernova picture, we 
suggest that the flow in PSl-lOjh underwent a shock or 
shocks at large radius from the BH, similar to a mech¬ 
anism prop osed for_ making superluminous supernovae so 
bright ('e.g.. lGal-Yanj|2012h ; see 33 In superluminous super¬ 
novae, it is believed that the high-velocity (~ 10, 000 km s^ 1 ) 
ejecta encounter a much slower wind expelled by the star 
earlier. In PSl-lOjh, perhaps a weak s hock at apocenter 
(IShiokawa et al.l 120151 : IPiran et al.l 120151 ') provided gas for 
the shell to shock on (see ( 18.2.21) . However, the slow speed 
of PSl-lOjh’s expanding shell (~ 1000 kms -1 ) would likely 
have hindered the effectiveness of such an external shock. 
We speculate instead that there could have been inter- 
nal shocks (similar to models for gamma-ray bursts; e.g., 
iRees fe Meszaro3ll994T ): deep inside the shell, perhaps more 
gas was accelerated outward from small radii (perhaps by 
radiation pressure), at speeds that (for some reason) were 
faster than the shell’s expansion. As this faster inner gas 
caught up to and collided with the slower-moving shell, the 
resulting shocks would convert the wind’s kinetic energy to 
radiation energy, producing the relatively large quantity of 
radiation observed (.Ek <C E l ra d). It would also impart mo¬ 
mentum to the shell, which could help explain why the ob¬ 
served expansion velocity was constant, despite the fact that 
the central black hole’s gravity should have continuously de¬ 
celerated the gas (see Figure [3]). This continual reheating of 
the gas could perhaps have maintained the temperature of 
the photosphere as well. Investigating these ideas further 
likely will require radiation hydrodynamical simulations. 


7 COMPARISON WITH SUPERNOVAE 

iGezari et al.l (j2012l ) argue against a supernova interpreta¬ 
tion for PSl-lOjh because of its long-lasting UV emission: 
GALEX still observed UV emission fully 375 (rest-frame) 
days after the light curve’s peak. They argue that a core¬ 
collapse supernova would have cooled to ~ 6000 K a month 
after the explosion; the lack of recent star formation in the 
host galaxy also makes a core-collapse supernova unlikely. 
Here we offer a further comparison between the properties 
of PSl-lOjh and those of supernova (SN) explosions. 


7.1 Photosphere and energy 

In ( 12.11 we calculated the radius of the photosphere of PSl- 
lOjh by comparing the temperature with the blackbody lu¬ 
minosity (eq. [3); we found that R p h increased with time 
for ~ 67 days, reached a maximum of R p h ~ 6 x 10 14 cm, 
and then notably decreased with time (Figure [lj. Because 
the effective temperature remained roughly constant, the 
expanding and receding photosphere caused the light curve 
to rise and fall. This scenario differs strongly from super¬ 
novae, where the radius of the photosphere generally in¬ 
creases with time up through the peak of the light curve 
and beyond; supernova photospheres only begin to recede 
once the whole flow has be come optically thin (e.g., [Arnett] 
1 19821 : IChornock et alJl20l4 . and references therein). For su¬ 
pernovae, the fading light curve is not due to a reced- 
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ing photosphere, but rather to the temperature of the gas 
falling ( mostly b ecaus e of adiabatic expansion of the gas). 
IChornock et al.l (12014 1 calculate R p h for PSl-llaf and find 
that its photosphere rises and falls as well; they also point 
out that such behavior clearly distinguishes that event from 
supernovae. 

This difference presumably comes about because the 
outflowing gas in PSl-lOjh is less massive (~ O.OIMq; eq. 
m Tan the e jecta in a typical supernova explosion (j> IMq; 
Smarttll2009l l: the column density in the expanding shell falls 
enough to become optically thin much sooner in PSl-lOjh 
than in a SN. Another difference is that in a supernova, the 
photosphere begins to recede at late times when it has cooled 
enough that the gas recombines and Thomson scattering 
ceases to provide opacity. PSl-lOjh is not seen to cool, and 
so Thomson scattering remains important. 

We note also that while the peak luminosity and total 
energy emitted from superluminous supernovae (SLSNe: see 
next section) are similar to PSl-lOjh, the photosphere for 
PSl-lOjh is a factor of ~ 10 smaller than that for SLSNe 
(e.g., PTF09cn d has R p h ~5x 10 15 cm at peak luminosity; 
IQuimbv et al.l201ll) . This observational difference is because 
PSl-lOjh is significantly hotter than SLSNe, and so needs a 
smaller photosphere to produce the same luminosity and 
energy output. 

Also at odds with supernovae is our finding that the 
kinetic energy in PSl-lOjh is much smaller than the radia¬ 
tion energy ( Ek E^a: ^3.21 ) ; supernovae t ypically exhibit 
Arad Ek (e.g.. lArnettlll982 : ISmarttll2009l) . We suggested 
in 1 )6.21 that this difference arises because the radiation in 
PSl-lOjh was mostly converted to potential energy instead 
of kinetic energy, since photons had to do work against the 
gravitational potential of a massive BH (which they don’t 
have to do in a supernova). We also suggested in 1 )6.2.1 1 that 
PSl-lOjh may not e xpand adiab atically and homologously 
as SNe typically do dArnettlll982h ; kinetic energy may have 
been converted back to radiation energy through shocks in 
the expanding gas. This would also have kept PSl-lOjh from 
cooling as SNe do. 

7.2 Light curves, spectra, and super-luminous 
supernovae 

We compare the light curve of PSl-lOjh with those of super¬ 
novae in Figure [5] We focus particularly on a growing class 
of supernovae called “super-luminous supernovae” (SLSNe): 
similar to PSl-lOjh, these events reach peak luminosities 
~ 10 44 ergs -1 and radiate total energies ~ 10 51 — 10 52 erg, 
10 or more times brighter than “common” SNe (IGal-Yaml 
2012lh Data for PSl-lOjh is plotted in orange diamonds over 
Gal-Yaml (l2012l Vs Figure 1. 

_Figure [5] and comparison with spectra in lGal-Ya3 

(120121) show that PSl-lOjh differs substantially from SNe, in¬ 
cluding SLSNe. Despite its large bolometric luminosity that 
would classify it as a SLSN, PSl-lOjh falls below the thresh¬ 
old when only r-band light is considered, because PSl-lOjh 
is significantly hotter and most of its radiation is emitted 
at ultraviolet wavelengths; UV emission in SNe is typically 
suppressed by line blanketing as well. PSl-lOjh fades more 
slowly than most supernovae, looking most similar in de¬ 
cline to I Gal-Yaml l|2012h ’s SLSN-R class, which is believed 
to be powered by radioactive decay of nickel-56; however, 



Figure 5. Figure 1 from IGal-Yaml d2012j ). with data for PSl- 
lOjh overplotted (orange diamonds). All supernova data is in 
the observed R band (see [GaTYam] |2012 |) ; da ta for PS-lOjh is 
in the observed r band (’from lGezarTer a].11201d) . The dashed line 
shows the threshold for classification as a super-luminous super¬ 
nova (SLSN). PSl-lOjh appears fainter than SLSNe because it 
is hotter, so less of its radiation emerges in the r-band, despite 
a peak bolometric luminosity similar to SLSNe. PSl-lOjh fades 
more slowly than most SNe, looking most similar in light-curve 
shape to SLSN-R events (but these have completely different spec¬ 
tral properties; see text). 

the optical spectra of SLSN-R events show many metal ab¬ 
sorption lines, completely absent in the spectra of PSl-lOjh. 
SLSN-II events fade more quickly than PSl-lOjh, and have 
spectra that contain strong hydrogen emission lines (unlike 
PSl-lOjh); these events are thought to be explosions that 
take place inside thick hydrogen envelopes (either inside pre¬ 
vious expulsions from a dying star, or inside an incredibly 
bloated star itself) which are heated and ultimately allow 
the radiation to diffuse out. SLSN-I events, driven perhaps 
by a combination of radioactive decay and interaction with a 
thick (hydrogen-free) envelope, also fade more quickly than 
PSl-lOjh. They are typically very blue and luminous in the 
ultraviolet (like PSl-lOjh), but their UV emission lasts a few 
months, not a year (as PSl-lOjh does). Their spectra show 
no hydrogen, but they also show no prominent helium close 
to the light curve’s peak, and do show broad absorption lines 
due to metals, which PSl-lOjh lacks. 

In this framework, the best bet for a supernova inter¬ 
pretation would probably be an event similar to a SLSN- 
II, in which an explosion took place inside a thick shell of 
gas. However, the shell would have to be helium-rich rather 
than hydrogen-rich—perhaps a star that shed its hydrogen 
envelope longer in the past, and shed its helium envelope 
shortly before explosion. This is similar to, e.g. , the s cenario 
proposed for a set of SLSNe in iQuimbv et all (120111 ): how¬ 
ever, in their case, the spectra showed neither hydrogen nor 
helium (and instead showed other element J 4 ^!). Moreover, 

16 In the optical. IQuimbv et al.l 1201 il l’s SLSNe show strong O 
II absorption lines; PSl-lOjh was too hot for 0 ++ to be present 
in significant quantities, so its spectra show no O II. 
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the requisite shell proposed bv lQuimbv et all (|201 il l weighs 
many solar masses and is produced by the violent pulsations 
of an extremely massive dying star (M* ~ IOOMq), highly 
unlike the puny 0.01 Mq shell we infer for PSl-lOjh. 

Velocities inferred for SLSNe’s expanding shells and 
shockwaves through the shells don’ t obviously ma tch t he 
properties of PSl-lO.jh either: e.g., Quimbv et al.l d201 ill ’s 
events (type SLSN-I in lGal-Ya m 2012’s nomenclature) have 
photosphere expansion speeds of ~ 10 4 kms _1 rather than 
10 3 kms _1 , and SLSN-II events show Ha lin ewidths of ~ 
4000kms _1 (e.g., ISmith et al.ll2008l : lGal-Yamll 2012 l l , some- 
what slower than ~ 9000kms _1 observed for PSl-lOjh. We 
note that the ratio of velocities of PSl-lOjh’s He II 4686 line 
(likely the velocity of a line-driven wind) to its photospheric 
expansion is ~ 5, while this ratio for SLSNe is ~ 0.3. This 
may be because SLSNe quickly become too cool to continue 
driving a wind. 

For a final comparison, even the spectrum of an “un¬ 
usual helium-rich” supernova, belie ved to have lost its hy - 
drogen envelope prior to explosion dKawabata et al]|201Cl ). 
looks nothing like PSl-lOjh: the supernova shows He I lines 
rather than He II, and they are in absorption rather than in 
emission. 

Thus we conclude that PSl-lOjh differs significantly in 
multiple important ways from known types of supernovae. 


8 DISCUSSION 
8.1 Summary 

PSl-lOjh is a transient discovered by the Pan-STARRS 
Medium Deep Survey in 2010 and identified as a candi¬ 
date tidal disruption event. This event has one of the best- 
sampled light curves of a TDE to date, including the rise to 
peak luminosity followed over two months. Optical and UV 
photometry show that the temperature of the photosphere 
was ~ 3 x 10 4 K and remained fairly constant over time. Op¬ 
tical spectroscopy showed two He II recombination lines in 
emission, but no hydrogen (or any other) lines. We have in¬ 
vestigated the nature of PSl-lOjh, making arguments based 
on both the continuum and spectroscopic observations. Here 
we summarize and discuss our conclusions, and point out di¬ 
rections for future work. 

(i) The observed constant temperature T e ff ~ 3 x 10 4 K 
and the observed light curve, when combined, yield the ra¬ 
dius of the photosphere R v h as a function of time, assuming 
that the emission has a blackbody shape. The photosphere 
is seen to grow, reach a maximum of ~ 6 x 10 14 cm, then de¬ 
cline. The growth is close to linear in time, increasing with 
a velocity vo ~ lOOOkms -1 . 

(ii) We interpret the growth of the photosphere to indi¬ 
cate the expanding edge of a shell of gas. The photosphere 
reaches its maximum size when the shell has an optical depth 
of unity; thereafter, the shell is optically thin and we see 
through it to gas inside. This constraint allows us to esti¬ 
mate the mass of gas in the shell, ~ O.OIMq. 

(iii) The mass and velocity estimates provide an esti- 
mate for the kinetic energy of the gas, Uk ~ 10 47 ergs _1 . 
iGezari et al.l ()2012ih estimated the total energy emitted in 
radiation -E ra d,boi,out ~2x 10 51 erg. We showed that the lu¬ 
minosity in the observed bands (with no assumption about 


the shape of the spectral energy distribution, i.e., without 
assuming the emission is thermal) was ^ 6 x 10 5 ° erg. The 
ratio of kinetic to radiation energy is thus ~ 10 -4 . This is a 
factor of a million smaller than the typical ratio for super¬ 
novae, ~ 10 2 H mar H (1200911 . 

(iv) He II emission lines are produced by recombination 
of He ++ ions with electrons. Photoionization calculations 
demonstrate that the observed lower limit of the ratio of He 
II 4686 A to Ha implies that the emitting gas is helium-rich; 
the gas does not have solar composition. 

(v) Based on these photoionization results, a likely possi¬ 
bility is that the disrupted object was a red giant core whose 
envelope was previously stripped. Such an object is dense 
enough that it requires a small BH (Mbh <) 2 x 10 5 Mq) to 
disrupt the object outside the BH’s innermost circular orbit. 

(vi) The (half) linewidth of the He II line and contem¬ 
poraneous measurement of the radius of the photosphere 
together imply the mass of the central object, ~3x 10 5 Mq. 
This estimate is close to the BH mass we inferred by inde¬ 
pendent argument in the previous point. 

(vii) From the luminosity of the He II 4686 line, we infer 
that the emission measure of the gas producing the line is 
{EM) ~ 3 x 10 65 cm -3 , similar to the independent estimate 
of the emission measure we made based on the radius of the 
photosphere (from continuum observations). 

(viii) We infer that gas from the disrupted star fell back 
to the BH at a rate much larger than the Eddington rate for 
an extended period (hours to a day) following disruption. 
Photons were trapped in the highly optically thick gas, and 
radiation pressure drove the expansion of the gas shell; the 
shell was gravitationally unbound from the BH, but (sur¬ 
prisingly) only barely. Trapping of radiation ensures that 
the light curve of the event is far from simply proportional 
to the mass fallback rate (Figure [3: future investigations of 
tidal disruption must account for radiation pressure. 

(ix) Outside the photosphere to Thomson scattering, res¬ 
onance lines (such as He La) continued to trap a fraction of 
radiation and drive a further wind, slightly faster than the 
expansion of the shell. We encourage observers to seek ultra¬ 
violet spectroscopy of future tidal disruption candidates, and 
predict that their ultraviolet spectra will show broad reso¬ 
nance lines, blue-shifted and seen in absorption (or perhaps 
having a P-Cygni profile), as wo uld characterize an outward - 
moving wind (also predicted in lStrubbe fe QuataertllilOlll 'l. 
We also encourage contemporaneous X-ray observations, to 
better constrain the incident SED for modeling photoioniza¬ 
tion in the line-emitting gas. 

(x) The emitted radiation -E ra d,boi,out ~ 2 x 10 51 erg is 
only 10 times less than the energy released in the dis¬ 
ruption of a helium core, ~ 3 x 10 52 erg. This implies 
that the radiation we observe did not start out at a small 
radius ~ RT.core ~ 2 x 10 11 cm and then adiabatically 
cool as the gas expanded. Instead, the photons must have 
been generated close to the radius where we observe them, 
R p h ;> (few) x 10 13 cm. We suspect that the expanding shell 
shocked on other gas, in order to make up for the adiabatic 
cooling as gas expanded from small radii, and to produce 
such a large amount of radiation relative to kinetic energy. 
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8.2 Implications & questions for future work 

8.2.1 Energetics 

IStrubbe fe Quataertl d2009l l argued that many TDEs would 
lead to super-Eddington rates of BH feeding, and thus 
radiation-pressure-driven outflows would be likely. They 
considered a case they called “edge-dominated,” in which 
the photosphere would be trapped in the sharp edge of the 
flow and would initially expand with time, rather than re¬ 
cede with time through a broadly-distributed wind (see Fig¬ 
ure [2] above). The expanding shell we infer for PSl-lOjh is 
sim ilar to t his p icture, b u t with two significant differences: 
(1) IStrubbe fc Quataertl l2009li assumed that the outflow 
velocity would be ~ 0.1c and so the region where photons 
would not be trapped (the diffusion time is less th an th e ad- 
vectio n time) could be neglected; (2) IStrubbe fe Quataertl 
(120091) did not expect the temperature to remain constant, 
although they did find only a weak time dependence for tem¬ 
perature Teff oc t~ 7 / 36 , which could still be consistent with 
the data. 

Our model of PSl-lOjh as an expanding shell of gas 
leaves and raises several questions: 

(i) Why was the outflow so slow, having such tiny kinetic 
energy, only barely unbound from the BH f* 13.2l * 16.21) ? Was 
the physics of the shock and accretion somehow regulated 
to distribute the radiation energy in such a way that the 
maximum amount of gas would be able to escape? Could it 
be that the gas was given a distribution of velocities, but the 
fastest shells spread out so much that they quickly became 
optically thin (so invisible), and the photosphere always lay 
in the slowest (barely unbound) shell? 

(ii) Why was the temperature close to constant, despite 
the fact that the photospheric radius ranged by an order 
of ma g nitud e over the period of observation (Figure [TJ? 
iMilleil d2015l) proposes that constant temperatures may be 
explained by disk winds: more super-Eddington feeding rates 
produce stronger disk winds, so that the accretion rate onto 
the BH is only weakly dependent on the feeding rate. Could 
this idea explain PSl-lOjh’s temperature in the context of 
the other physical properties of the event we have inferred? 

(iii) Was the expansion of the shell non-adiabatic, as we 
argued in *16. 2. IP In that section, we suggested that there 
may have been an additional shock, converting kinetic en¬ 
ergy back to radiation. What caused this shock? Could this 
process also explain the observation that the gas was barely 
unbound, and the roughly constant temperature? 

We discuss these points further in the following subsec¬ 
tion. 


8.2.2 Orbit circularization & shocks at large radius 

IShiokawa et al.l ((20151 1 perform general relativistic hydrody¬ 
namic simulations of the disru ption of a white dwarf by a 
500Mq BH, followed up by IPiran et alj (|2015l l. These au¬ 
thors find that the process of circularization of orbits of 
stellar gas particles following disruption is slower and less ef¬ 
ficient than previously thought: the shock at pericenter may 
be weak, and a second and third shock close to apocenter of 
the most tightly-bound gas may be partially responsible for 
circularization as well. They find that the time for (partial) 
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circularization of orbits is ~ 3 — 10 times longer than pre¬ 
vious expectations (~ (3 — 10)tf a iiback: see eq. [52J, and the 
maximum feeding rate to the BH is thus ~ 10 times less. 
Although these effects would stretch the black curve in our 
Figure|5]for PSl-lOjh out in time, they do not stretch it out 
enough to make it match the green curve: the timescales for 
Mfaiiback and observed luminosity differ by a factor of ~ 10 3 , 
not just a factor of ~ 10. 

These authors indicate that the pericenter shock is 
strong if the disruption radius is close to the BH’s event hori¬ 
zon (i.e., if Rt <j 5R S ), since relativistic apsidal precession 
and potentially Lense-Thirring preces sion cause significant 
crossi ng of orbits (e.g., recent work by Bonnerot_et_aL 2015J; 


Ha yasaki, Stone fe Loebll2015l : iGuillochon fe Ramirez-Ruid 
2015ll . Provided Mrh > 10 5 Mm. the helium core disruption 


of PSl-lOjh falls in the regime where the pericenter shock is 
stro ng enough to make circularization rapid. 

IPiran et ah IHqH) suggest that the radiation seen in a 
tidal disruption event comes from the shocks near apocenter 
rather than the shock at pericenter (or from an accretion 
disk at yet smaller radii). However, we find that this does 
not work for PSl-lOjh: the observed low mass (~ O.OIMq) 
combined with the low orbital velocity at i? ap0 (see eq. 1211 ) 
shows that the energy liberated by the apocenter shocks is 


2GM BH M she ii 

R&po 

2 X 10 50 ( ^BH \ f Afshell / flapo \ 

^ 2 x 10 5 Mq ) O.OIMq J V 6 x 10 12 cm) 


an order of magnitude too small to explain the observed 
radiation energy. Nevertheless, we found in * )6.2.1l that the 
expanding shell likely experienced a shock at large radius 
that produced most of the observed radiation: we offer the 
speculation that perhaps gas that circularized through the 
apocenter shock was partially responsible (see * 16.2.111 . The 
process of circularizing the orbits of disrupted stellar ma¬ 
terial is clearly complicated and may be important for un¬ 
derstanding PSl-lOjh and other disruption candidates. We 
defer further investigation of these ideas to future work. 


8.2.3 Photoionization calculations 

We reiterate that if ionized hydrogen had been present in 
significant quantities, it would have recombined in steady 
state to produce strong Balmer emission lines, which were 
not seen. Recombination times were short (eq. lB2l) ; every 
hydrogen photoionization would have been balanced by a 
recombination cascade. Arguments in the literature that 
“over-ionization” of hydrogen provides a means of suppress¬ 
ing Balmer emission do not make sense. The only way to 
suppress Balmer emission is through absorption by excited 
neutral hydrogen, but velocity gradients in the outflowing 
gas suggest that Ha would not have been optically thick, 
and so such absorption would not have been significant. We 
conclude that PSl-lOjh comprised gas that was helium-rich, 
not of solar composition. See Appendix 0 for detailed ex¬ 
planation. 


(34) 

erg, 
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8.2.4 Mechanism to produce helium core 


We argue in Appendix [A] that the spectrum of PSl-lOjh 
is inconsistent with the disruption of a main-sequence star, 
and argue instead for the disruption of the helium core of a 
stripped red giant star. How could this come about? Detailed 
investigation of this crucial question is beyond the scope of 
this work, but we briefly offer a few comments here. 

First, could tidal stripping of a red giant 
star leave behind an object whose disruption 
cou ld produce the requisite line ratio 1/4686/Lhq? 
iMacLeod. Guillochon fe Ramirez-Ruid (120121 ) study the 
tidal stripping by a BH of mass Mbh = 10 6 Mq of a 1.4Mg 
star at various stages of post-main-sequence evolution. 
Their RG I model has a core of mass M core = 0.28 Mq and 
envelope of mass M env = I.IMq. The maximal stripping of 
the envelope (by the deepest disruption) removes ~ 70% 
of the envelope mass, leaving ~ 0.3Mq still gravitationally 
bound to the core. If we approximate that the envelope 
is purely hydrogen and the core is purely helium, and 
that the species fully mix following disruption, we find 
that this scenario leads to a hydrogen mass fraction 
X = Mh/(Mh + Mii e ) ~ 0.5 (c.f. eq. IA3I) . or a number 
density ratio 


riHe Mfje 
Hh 4Mh 


(35) 


Equation (lA2l) shows that in the optically thin limit, 
T4686/Thck oc nne / nn . Our Cloudy results in Figures lATI 
and IA2I (top left parts of Panels (a) and (d)), which as¬ 
sume a density ratio nHe/nu =0.1, show that the line ratio 
under optically thin condition£3 is ~ 0.4 — 1.8. This indi- 
cates that to reach the observed line ratio I/ 4686 /Tho > 5 
dGezari et al ] 2013), the number density ratio would need to 
be ~ 3 — 10 times as high as assumed by the Cloudy calcula- 
tio ns, or nn e /nn ~ 0.3 — 1. (If we make the same argument 
for lGaskell fe Roias Lobosll2014l ’s result of L 4686 /Ehci > 3.7, 
we find )ih 8 /)ih ~ 0.2 — 0.9.) On the low end of these re¬ 
sults (0.2 —0.3, which we derived for the incident SEDs clos¬ 
est to that of PSl- lOjh), the number density ratios come 
close to those of IMacLeod. Guillochon fc Ramirez-Ruid 
d2012h’s stripped red giants (eq . l35l) . This suggests 
that IMacLeod! Guillochon fe Ramirez-Ruid d2012l Vs strip¬ 
ping mechanism could be a viable one for producing a 
helium-rich core whose disruption spectrum would look like 
PSl-lOjh’s; complete stripping of the hy drogen envelope 
may not be required. Calculations bv lCheng fc Evansl (l2013i j 
also indicate that the outer layers of the (mostly) stripped 
core are tidally heated each time the core passes pericen- 
ter around the BH, and thus the remainder of the stellar 
envelope could be removed after several passes. 

What happens a fter the envelope of the red giant is 
tidally stripped away? IBogdanovic. Cheng fc Amaro-Seoanel 
d2014h propose the following scenario. Following the tidal 
stripping, half of the envelope is accreted and half es¬ 
capes unbound, as would happen in full disruption of a 
solar-type star, though with an abrupt cut-off in accre- 


17 The higher values than in eq. lA2l l are due to the temperature 
dependence of the recombination coefficients. We discuss the op¬ 
tically thin limit because velocity gradients in the outflowing gas 
likely reduce optical depths dramatically; see j|A5l 


tion at late times b ecause the core gravitationally retains 
the g a s closest to it (MacLeod. Guilloc hon fe Ramirez-Ruid 
12012 : iGuillochon fc Ramirez-Ruid I2013T ). (This accretion 
may have produced a bright flare, but evidently not dur¬ 
ing any period when this region of the sky was being ob¬ 
served.) The stripped core then gradually becomes tidally 
heated by the BH’s gravitational field, while its orbit 
slowly circularizes and decays by the emission of grav¬ 
itational waves (as an EMRI: extreme mass-ratio inspi- 
ral) . For just the rig ht or bital a nd stel lar parameters, 
IBogdanovic, Cheng fe Amaro-Seoanel (j2014h find that the 
core will retain (rather than radiate away) its heat long 
enough to expand and “lift its degeneracy” enough to be 
tidally disrupted by the BH before gravitational waves carry 
it through the ISCO or event horizon. Note that this scenario 
considers a BH of mass Mbh i> 10 6 Mq, which is too massive 
to tidally disrupt a helium core that isn’t tidally heated; a 
lower-mass BH (Mbh $2x 10 s Mq; eq. 1191 can do the job 
without tidal heating. 

In addition to hypothesizing scenarios to tidally 
disrupt helium-rich objects, it is also critical to con¬ 
sider how common such scenarios are likely to be. 
IMacLeod. Guillochon fe Ramirez-Ruij |2oT3) make careful 
estimates of disruption rates of stars at different stages 
of evolution (see their Figure 13 in particular). Main- 
sequence stars are by far the most common disruptions, 
but red giant disruptions contribute ~ 10% to the to¬ 
tal, especially for stars slightly more massive than the 
sun and above. (This work did not consider BHs of low 
enough mass [Mbh <[ 2 x 10 5 Mq] to disrupt red gi¬ 
ant core s.) Future researc h cou l d com bi ne an d exte nd the 
wor k of IMacLeod, Guillochon fe Ramir ez-Ruij (120121 ) and 
IBogdanovic, Cheng fc Amaro-Seoanel ( 20141) to determine 
theoretically how often the full scenario takes place: a red 
giant’s core is tidally stripped of its hydrogen envelope, then 
it spirals in to the BH via gravitational wave emission, and 
then is ultimately disrupted. The theoretical ratio of this 
rate to the rate of all disruptions would be illuminating for 
deciding how plausible this scenario is for explaining PS1- 
lOjh. 


8.2.5 Comparison with other events 


As the number of tidal disruption candidates with de¬ 
tailed observations grows, preliminary comparisons of 
rates can become possible. As mentione d above, PTF09ge 
llKasliwal et al.l[2009l : lArcavi et al .1 EfoTih is a recently pub¬ 
lished candidate that also showed broad He II 4686 emis¬ 
sion and no hydrogen. Could this event also be the dis¬ 
ruption of a stripped helium core? Would the discovery of 
two such events be in line with theoretical rate expecta¬ 
tions? Other tidal disruption candidates with optical spec¬ 
troscopy show both heli um and hydrogen or emission, or 
only hydrogen emission (lArc avi et all 120141 : iHoloien et alJ 
1201-4 |yaJL(Velzen_et_ai] |201l|), or no opt i cal lines at all 
( Cenko et al. 2012i : Ivan Velzen et al~]|201ll : IChornock et al.l 
120141) . 


The type of analysis we have done for PSl-lOjh required 
observations of the rise of the light curve, ultraviolet pho¬ 
tometry, and spectroscopy, in order to derive the tempera¬ 
ture and thus the photospheric radius, and to constrain the 
composition of the disrupted object. We encourage transient 
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surveys to seek all of these data for candidate disruptions: 
it will be very useful to know if other events also show a 
clearly expanding then contracting photosphere, how fast it 
is moving, how massive is the gas, and the energetics of the 
event; and also to perform photoionization calculations to 
compare with spectra that can give clues to the composi¬ 
tion of the disrupted object (and thus perhaps its density 
and constraints on the BH mass). 


8.2.6 Low-mass black hole 

We estimated that the mass of the BH was Mbh ~ 2 x 
10 5 Mq or less, in order for the BH’s tidal forces to be able 
to disrupt a dense core of a red giant outside the BH’s ISCO. 
(We independently estimated that Mbh ~ 3 x 10 s Mq by 
comparing the (half) linewidth of He II 4686 with the radius 
of the photosphere.) Such modest BHs are expected to be 
very common in the universe, but are difficult to find via 
stellar dynamics or AGN emission, because the radius of 
influence and Eddington luminosities are so small. 

PSl-lOjh may offer an exciting firsl[3 significant BH 
mass constraint using a tidal disruption event. The host 
galaxy of PSl-lOjh has a total stellar mass around 3 x 10 9 Mq 
but unknown morphology; SDSS contains an image of the 
galaxy, but with too poor resolution to decompose a bulge 
from a disk if present—so the bulge could be ~ 10 times less 
massive than the total stellar mass. At these low masses, 
the Mbh — Mb u i ge relatio n increases in scatter and becomes 
poorly constrained (e.g., iGreene. Ho fe Barthll2008l f. A BH 
of mass 2 x 10 5 Mq in a host bulge of 3 x (10 s — 10 9 )Mq 
appe ars rou ghl y in line wit h pre vious Mbh — Mb u i ge result s 
(e.g., IGreene. Ho fc Bartbll2008l ; [j iang . Green e & Holl201l|j . 
PSl-lOjh could offer a new datapoint for the A/bh — Mb u i ge 
and M — a relations if additional observations (decomposing 
the bulge from the disk, and measuring the velocity disper¬ 
sion) of the host are made. 


8.2.7 Super-Eddington accretion physics 


Our BH mass estimate Mbh ~ 2 x 10 5 Mq led us to conclude 
that the rate of feeding the BH was much higher than the 
Eddington rate, by a factor of ~ 10 4 close to the peak (eq. 
m Figure [I]). 

Such highly super-Eddington feeding generates a large 
amount of radiation and is expected to produce a compli¬ 
cated flow, with photon trapping, photon diffusion, outflows, 
accret ion, and possibly o t her p r ocesses all playing roles 


accret ion, and possi bly ot her p r ocesses all playm; 
(e.g., Shakura fc Sunyae vl 19731: Abramo wicz et ahl 


iBegelmanll 19791 : lOhsuea fc Mineshigell2007l ~). Past studies of 
such flows have typically inferred low efficiencies for the con¬ 
version of mass feeding to escaping energy, because the gas 
becomes extremely optically thick and the diffusion of radia¬ 
tion out is ineffective; most of the radiation remains trapped 
and is soon swallowed by the BH along with the gas in which 
it is entrained. _ 

Ijiang, Stone fc Davisl ||2014| ) recently performed a 


18 Other events have claimed Mbh measurements based on their 
light curves, but we believe that current understanding of the 
precise physical mechanisms producing the light curves is still 
too preliminary to make such claimed measurements. 


global 3-dimensional radiation magneto-hydrodynamic 
(MHD) simulation of super-Eddington accretion and a find 
significantly higher efficiency r/esc = Et ot ,out /Me 2 = 0.045. 
Their simulation finds that the gas is turbulent, and cooling 
is dominated by magnetic buoyancy rather than photon dif¬ 
fusion; as a result, a significantly higher fraction of energy 
is able to escape. They also find that the escaping energy 
is divided roughly 5 : 1 between radiation and ki netic en¬ 
ergy: in our terminology, Ek ~ 0.2E r *d ^..t. lMcKinnev et ahl 
d2014f) recently performed a 3D general-relativistic radia¬ 
tion MHD simulation of super-Eddington accretion onto a 
rapidly rotating BH. They find that about 1% of matter 
fed to the BH is converted to radiation (?? ra d ~ 0.01) and 
about 80% of matter fed to the BH is actually accreted 
{Va.cc ~ 1 — 0.8 = 0.2 in their nomenclature). 

It is interesting to compare these simulation results with 
accretion properties we infer for PSl-lOjh. The overall ef¬ 
ficien cy we find f or P Sl-lOjh is rfesc ~ 0.02 ( N6.2I) . simi¬ 
lar to I Jiang. Stone fc Davisl ll2014l Vs result. Evidently, a re¬ 
markably high fraction of the mass of the disrupted star 
is converted to energy which manages to escape to great 
distances. However, we find that it does not escape imme¬ 
diately, but rather is trapped for 10 2 ^_10 3 _times the in¬ 
flow (fallback) time. Ijiang, Stone fc Davisl (l2014| j simulated 
a domain corresponding to size r -C _R ap o, so that their 
simulation could not account for any (potentially optically 
thick) post-apocenter shocked gas, as it is well outside the 
radius of their disk. We also infer a much lower kinetic en¬ 
ergy (Ek ~5x 10 -4 E ra d, O ut) than they do, but again, their 
simulation cannot address the effects of gas infalling from 
Rapo- They find that the mass-loss rate from the disk is 
~ 30% of the mass accretion rate onto the black hole, while 
we infer a ratio (on much larger scales) of ejected mass to 
accreted mass of ~ O.O lMf?) /0.15Mp> ~ 7%. 

Comparing with iMcKinnev et alj 120141) . we infer 
Viad ~ 2E ra d,boi,out/M core c 2 ~ 0.007, close to their value 
of 0.01 (though ours is integrated over time, rather than 
comparing at a specific moment). They find that their wind 
carries away some 20% of t he mass , midway between the 
result of ljiang. Stone fc Davisl (|2014| ) and our inferred ratio 
for PSl-lOjh. 

Future super-Eddington accretion simulations may fo¬ 
cus more on tidal disruptions by starting from initial con¬ 
ditions given by studies of orbital circularization of falling- 
back stellar debris (see references in 48.2.21) . Continued com¬ 
parisons between observed properties of tidal disruption can¬ 
didates and accretion simulations may give further insights 
into accretion physics relevant in other contexts as well, such 
as ultra-luminous X-ray sources and the growth of massive 
BHs at high redshift. 
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APPENDIX A: PHOTOIONIZATION 
CALCULATIONS: NOT THE DISRUPTION OF 
A MAIN-SEQUENCE STAR 

In this section we show that the observed spectrum of PS1- 
lOjh is unlikely to be produced by the disruption of a 
main-sequence star. We perform calculations using the pub- 
licl y availab l e photoion ization code Cloudy (last described 
by iFerland et al.l l2013l ~) , and compare calculated line ra¬ 
tios of He II 4686 to Ha with the lower limit observed for 
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PSl-10jh. iGezari et al.l (120121 ') find that any Ha emission is 
at least 5 times fainter than the He II 4686 A emissio 10 
(±4686 /± Hn > 5);IOask ell & Rojas Lobosl d20l4 ) separately 
analyze iGezari et al.l (I2012T) ’s spectrum and suggest that 
±4686/±Ha > 3.7 ±25%. 


A1 Review of photoionization physics 


We begin by briefly reviewing the physics of photoionization 
and how abundance limits may be inferred from observed 
line ratios. Photons incident from a central source (perhaps 
an accretion disk close to the BH’s event horizon) strike 
a cloud of gas some distance away from the source; if the 
photons are energetic enough, they photoionize atoms and 
ions in the cloud. After a recombination time t rec feq. IB2I) . 
electrons and ions begin colliding and recombining, generally 
to an excited state of the atom or ion; the electron then 
spontaneously falls to lower and lower energy levels, emitting 
line photons with each transition. In an optically thin region 
(from which all line photons can escape), the rate of emission 
of a given line transition is 

±iine ~ »i»e«!™(n, T)V (hv i ine ), (Al) 


where rii and n e are respectively the ion and electron num¬ 
ber densities, are" 6 is the effective recombination coefficient 
of the transition, V is the emitting volume, and hv ii ne is the 
energy of the transition. He II 4686 A and Ha are lines of 
singly ionized helium (4 —> 3) and neutral hydrogen (3 —> 2), 
with well-known recombination coefficients. For example, at 
T = 2 x 10 4 K, 04686 = 1-77 x 10~ 13 cm 3 s -1 is the effec¬ 
tive Case B recombination coefficient for He II 4686, and 
Oho = 6.04 x 10~ 14 cm 3 s^ 1 is th e effective Case B recom¬ 
bination coefficieni 23 ! for Ha (e.g.. [Osterbrockill989l ). (Note 
that the ratio of these coefficients is 2.9; more discussion of 
the temperature dependence of this ratio is below in j)A4l ) 
The line flux ratio of He II 4686 to Ha in an optically thin 
region is therefore approximately 

±4686 ^ n e n He ++ 04686 (hv 4686) 

Ifift n. e 71]4+ OHck (tllS H a ) 

~ Q 1 ^ n He++ / W H+ ^ ^ Q4686 /QHa ^ 


If the following conditions are met—(1) hydrogen and 
helium are both highly ionized (so that riH ~ n H + and nn e ~ 
n He ++), (2) the region is optically thin to Ha and 4686 A 
photons, and (3) the temperature of the gas is relatively low 
(T -C 10® K so that 04686 /oh« ~ 2.9)—then the line ratio 
offers a direct estimate of the abundance ratio of helium to 
hydrogen. The line ratio lower limit ±4686 /Iho > 5 implies 
riHe/riH >1.2, which corresponds to a mass fraction 


X = 


n H 

n H +4n He 


< 0 . 2 , 


(A3) 


19 IGezari et al.1 (|2012h published no error bars for this number, 
because of the difficulty in subtracting off the unknown bright¬ 
ness of the host galaxy. Now that the event has faded below the 
brightness of the host, it is possible to determine more precise 
constraints on L 4686 /^Ha > as Gezari et al. are currently working 
on (personal commun ication). 

20 lOsterbrockl i 19891) gives a h/ 3 = 1-62 x 10 14 cm 3 s 1 and 
jua/ja/3 = 2.76, both at T = 2 x 10 4 K, where j H x oc «hx"Hx- 


far b elow the ab un dance ratio X « 0.7 found in solar-type 
stars. ICezari et al make essentially this argument 24 ! 

and propose on this basis that PSl-lOjh is the disruption of 
a helium-rich object, such as the helium core of a stripped 
red giant star. 

These simplifying conditions are generally found for rel¬ 
atively high fluxes (to highly ionize the gas), low densities (to 
keep lines optically thin) and soft SEDs (to keep the gas rela¬ 
tively cool). When these simplifying conditions are not met, 
other additional effects become important, e.g., collisional 
ionization, 3-body recombination, and line absorption. Nu¬ 
merical calculations, such as those using Cloudy, become 
crucial for capturing these complicated radiative effects. 


A2 Proposal that PSl-lOjh is disruption of a 
main-sequence star 

ICuillochon. Manukian & R.amirez-R.uij ( 2014 1 propose that 
PSl-lOjh is the disruption of a main sequence star: they pro¬ 
pose that hydrogen is indeed present in cosmic abundances 
in the outflowing material, but radiative transfer effects sup¬ 
press emission from hydrogen relative to helium below de¬ 
tectab le limits. They appeal to Figure 1 in lKorista fe Goadl 
d2004i), w ho use Cloudy to calculate theoretical equivalent 
width^j (EWs) of various emission lines in the spectra of ac¬ 
tive galactic nuclei (AGN). The Cloudy calculations assume 
cosmic abundances of hydrogen and helium (riHe = OTtih), 
and are performed for different combinations of hydrogen 
gas density bh and hydrogen-ionizing number flux 


$ = 


i ru, 

47r r 2 I hv V ’ 

J XH 


(A4) 


with Xu = 13.6 eV the ionization potential of hydrogen. 
In the region of this figure where bh ~ 10 12 cnR 3 and 


a few A, while the Ha EW is < lA (although it is not clear 
from the figure that He II 4686 is fully 5 times brighter than 
Ha even in this region). Guillochon et al.’s idea is that PS1- 
lOjh was in this region of parameter space during the two 
epochs when spectra were observed: the y assert that the lack 
of observed Ha is consisten t withlKorista fe Goadl (120041 ') ’s 
calculated line ratios. Since iKorista fe Goadl ( 20041 ) ’s calcu¬ 
lations assume cosmic abundances of hydrogen and helium, 
as would be present in a main sequence star, Guillochon et 
al. assert that PSl-lOjh is the tidal disruption of a main 

sequence star. _ 

_However, iGuillochon. Manukian fe Ramirez-Ruij 

(l2014ll draw these conclusions without taking into account 
the magnitude of the observed lower limit of He II 4686 
to Ha emission (±4686 /±h« > 5), and without performing 
their own Cloudy photoi onizat ion c alcu l ations fo r th e 
relevant input parameters. [Gaskell fc Rojas' Lobosl d2014h 
do perform Cloudy calculations to investigate the line ratio 
for PSl-lOjh, and find that the ratio ± 4686 /±Ha can be 
~ 4 for a narrow range in density around ~ 10 44 cm~ 3 . 


21 In lGezari et al.l d2012h Supplementary Information §7; we have 
fixed the small mistake where He+ and H° were written instead 
of He++ and H+. 

22 These equivalent widths are relative to the continuum flux at 
1215 A. 
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They explain that this occurs when Ha becomes optically 
thick and the line ratio approaches its blackbody limit (see 
below our eq. lA5l and surrounding discussion). However, 
as explained further below, (1) these Cloudy calculations 
assume a hard AGN spectral energy distribution when 
a thermal SED is more appropriate, and (2) it is likely 
that these calculations are not converged. Previous work 
has also not considered how velocity gradients in the gas 
dramatically reduce line optical depths and likely reduce 
maximal line ratios, which is likely a serious problem for all 
of these results so far (see am 

We study these issues in detail, and find that there are 
significant problems with the interpretation that PSl-lOjh 
is the disruption of a main sequence star. 


A3 Cloudy calculations 


We have performed suites of photoionization calculations 
using the latest rele ase of Cloudy C13 (version C13.O10) 
jFerland et al.l 120131 ') . looking for regions of parameter 
space where the He II 4686 line luminosity is more than 
5 times the Ha line luminosity. Note that Figure 1 in 
iKorista fe Goad! (120041 ') was created from calculations us¬ 
ing significantly older Cloudy version v90.04. The up¬ 
dated version of Cloudy includes a much more sophisti¬ 
cated treatment of the hydrogenic and helium-like atomic 
energy levels, giving more accurate line luminosities for 
transitions between these energy levels (K. Korista, per¬ 
sonal communication), crucial for investigating Z/ 4686 /£hc« 
in PSl-lOjh. Unfortunately, however, in versions to date, 
Cloudy has difficulty converging as the number of “resolved” 
atomic energy levels increases, in the region of parame¬ 
ter space where both the density and incident flux are 
high (R. Po rter, personal communication). Althou gh not 
disc ussed in [Guillochon, Manuk ian fe Ramirez-Ruid (j2014h 


and iGaskell fe Rojas Lobosi (|2014h . these issues should 
present problems for conclusions in these papers. We present 
our Cloudy results with caveats about convergence where 
relevant. 

An important consideration in Cloudy calculations is 
the spectral energy distribution (SED) of the source whose 
light is incident on the theo retical cloud of gas who se emis¬ 
sion lines we wish to model. IKorista fc Goadl (|2004l Vs calcu¬ 
lations are all for a hard active galactic nucleus SED includ¬ 
ing a luminous X-ray taiQ- Guillochon et a l.’s assertion that 
PSl-l Ojh had solar composition rests on IKorista &; Goadl 
(l2004h ’s calculations using this incident SED, but the ob¬ 
served SED in PSl-lOjh is quite different from that of an 
AGN. PSl-lOjh was followed up by Chandra X-ray Observa¬ 
tory 315 days after peak, and showed an upper limit of only 
Lx(0.2 — 10) keV < 5.8 x 10 41 ergs -1 . Incident X-ray power 
is important because X-rays heat the g as, which leads to 
reduced recombination coefficients (e.g., IVerner fe Ferlandi 


23 Most recent Cloudy version C13.03 has no revisions strongly 
relevant to our calculations. 

24 The command in Cloudy is “agn 6.683, -1.2, -1.2, -0.9” (K. 
Korista, personal communication), where the numbers in order 
indicate the logarithm of the temperature of the “Big Bump” 
component, the X-ray to UV ratio o 0 x, the low-energy slope of the 
Big Bump continuum cqjV> an d the slope of the X-ray component 

o X - 


Il996h and also smaller optical depths in lines. Gezari et al. 
fit the observed SED with single blackbody curves at effec¬ 
tive temperatures of Tefj^ow = 2.9 x 10 4 K and Left,high = 
5.5 x 10 4 K (depending on the unknown level of extinction). 

We perform Cloudy calculations for four incident SED 
shapes, and a range of radiation fluxes and densities. For 
each model, we perform calculations using n res = 10, 15, 25, 
and 35 “resolved” atomic energy levels for hydrogen and he¬ 
lium, to check convergence: when a calculation is converged, 
increasing the number of resolved energy levels should not 
change the results. We calculate the ratio of output fluxes in 
an individual line between the n les = 25 and 15 calculations 
(“Ratio I” for Ha and He4686), and between the n res = 35 
and 25 calculations (“Ratio II” for Ha and He4686). In some 
cases (especially at high density and low incident flux), the 
calculations for n les = 35, and occasionally for n les = 25, did 
not complete in 48 hours (on an Intel Xeon 1.6 GHz proces¬ 
sor). We adopt as our criterion of convergence that (1) at 
least nres = 25 must have completed successfully, and (2) if 
Tires = 35 completed successfully, then Ratio II must lie be¬ 
tween 0.9 and 1.1 for both Ha and He4686; if only n res = 25 
completed, then Ratio I must lie between 0.9 and 1.1 for 
both Ha and He4686. Otherwise, we consider the calcula¬ 
tion not to have converged. For converged calculations, we 
report values from n res = 35 if available, or from n re s = 25 
if not. _ _ 

We first recalculate IKorista fe Goadl d2004h ’s Figure 
1 (on which Guillochon et al. base their claim of so¬ 
lar composition). We perform calculations on a 10 X 8 
logarithmically-distributed grid of density and incident flux, 
with n h G IIP 7 ,10 14 j cm -3 and $ G [10 17 ,10 26 ] s -1 cm -2 ; as 
IKorista fe Goadl (|2004l ) do, we assume cosmic abundances 
and a column density of Ah = 10 23 cm -2 . When highly ion¬ 
ized, this column will have an optical depth to Thomson 
scattering r es ~ 0.1, safely below the order-unity limit for 
which Cloudy is designed. We present our results in Panel 
(a) of Figure lAll Each colored square^ represents the line 
ratio L 4686 /Lhc calculated with Cloudy; white squares rep¬ 
resent regions of parameter space where the calculations did 
not converge (as described above). For this AGN incident 
SED, we find that the highest value of L 4686 /Lh<* is only 
2.8 (found for n H = 10 12 cm -3 , 4> = 10 23 s -1 cm -2 ). 

Next we perform Cloudy calculations for thermal SEDs 
at T = 2 x 10 5 K, 5.5 x 10 4 K and 3 x 10 4 K, at the same 
range of radiation fluxes and densities as the AGN calcula¬ 
tion above. We present our results in Figure 1X11 Panel (d), 
and Figure [A~2l Panels (a) and (d), and find that the high¬ 
est value of L 4686 /Lhck is 3.4 (found for T = 2 x 10 s K, 
nH = 10 12 cm -3 , $ = 10 24 s -1 cm -2 ). 

All of these results suggest that large densities and large 
ionizing fluxes are the most likely conditions to produce sig¬ 
nificant line ratios L^se/Lna (as explained more in i )A4l) . 
but also that the observed line ratio lower limit of 5 is too 
large to have been achieved with cosmic abundances for ei¬ 
ther an AGN or a thermal SED (of appropriate tempera¬ 
ture) . 


25 Note that Figures lAll and [A2l appear with different color tables 
in the published MNRAS version, at the request of the editor 
for better black-and-white printing. The data represented in the 
figures here and those in the published version are the same. 
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Observations of the PSl-lOjh continuum can indicate 
the relevant region of ( n , <3>) parameter space at the time 
of the first spectrum. As described in m\ and Footnote 
[5] if we assume that the gas has cosmic abundances, we 
infer that the number density in the shell at that time was 
respect ^ 1 x 10 cm for Xeff ■ 5.5 x 10 K, or rispect ^ 
8 x 10 10 cm -3 for T e g = 3 x 10 4 K. For an incident SED 
with temperature T e g = 5.5 x 10 4 K, the hydrogen-ionizing 
number flux is $ = 1 x 10 25 s -1 cm -2 , and for T e g = 3 x 
10 4 K, the flux is <f> = 4 x 10 23 s -1 cm -2 . Based on these 
estimates, the regiorFH in the suites of Cloudy calculations 
most relevant to the observations of PSl-lOjh is thus nu ~ 
10 11 cm -3 and $ ~ 10 25 s -1 cm -2 . For these parameters, 
our Cloudy calculations (which are all converged) find line 
ratios far below the observed lower limit of 5: 1/4686 /Ihh = 
1.1,1.5, 0.6,1.2 for thermal SEDs with T = 3 x 10 4 K, 5.5 x 
10 4 K, 2 x 10 s K, and the AGN SED. 

A4 Interpretation of Cloudy results 

Understanding the physical processes important in the dif¬ 
ferent Cloudy calculations will help us decide how to apply 
them to understanding PSl-lOjh. We begin by explaining 
the results for a thermal incident SED at T = 2 x 10 s K de¬ 
picted in Panels (d) - (f) in Figure [All and then explain how 
the results differ for different incident SEDs. We examine the 
different regions of ( n , 4>) in turn. 

The top left corner of these plots, where the density 
is low and the incident ionizing flux is high, corresponds 
to the simple conditions described around equation (IA2I) 
in m the large ionizing flux ensures that hydrogen and 
helium both are almost completely ionized throughout the 
region; meanwhile, the low density keeps the recombination 
rate low, so the populations of neutral hydrogen and singly 
ionized helium are tiny, making the region optically thin to 
both the Ha and He II 4686 lines (Panels (e) and (f)). The 
line ratio (Panel (d)) is thus close to the simple theoretical 
value of 0.4 described in equation (IA2I) . 

In the bottom part of these plots, the ionizing flux is 
too small to ionize most of the cloud 27 !. The majority of 
hydrogen and helium is neutral, with little H + around to 
recombine and produce Ha, but even less He ++ around to 
recombine and produce He II 4686 emission, since the ion¬ 
ization potential for helium is higher than for hydrogen (and 
helium requires two ionizations to become He ++ rather than 
one, as hydrogen does to become H + ). This explains the tiny 
value for 1/4686/ Luc, in the bottom part of the figure. (The 
bottom part of Panels (e) and (f) shows that Ha is optically 
thick but He II 4686 is optically thin— because He II absorp¬ 
tion would require a substantial population of He + ions, but 
almost all of the helium is neutral.) 

In the top right part of the plots, the ionizing flux is 
large enough to keep the region almost completely ionized, 

26 Note that iGuillochon, Manukian fe Ramirez-Ruijj ll2014l 'l 
identify a somewhat similar region of parameter space for the 
time of the first spect rum, tin ~ 10 12 — 10 13 cm~ 3 and 4> ~ 
10 23 — 10 24 s -1 cm~ 2 . [Caskell fe Rojas Lobosi ll2014h also iden¬ 
tify riH ~ 10 11 cm 3 as leading to the largest line ratios, although 
they do not specify the corresponding flux. 

27 This situation is sometimes called “radiation bounded;” its 
opposite is called “matter bounded.” 


yet the density is large enough that a relatively high recom¬ 
bination rate maintains small but significant populations of 
neutral hydrogen and singly ionized helium. Under these 
conditions, the Ha and 4686 A lines can become optically 
thick, and processes such as collisions and 3-body recom¬ 
bination become important. The limiting behavior, as the 
cloud becomes optically thick at all wavelengths, is that the 
line ratio tends towards the ratio of blackbody flux at the 
line wavelengths, i.e., 


F/4686 ^ f 4686 A \ _ g 

L Ha ~ \6563~A ) 


(A5) 


if these wavelengths are on the Rayleigh-Jeans tail (where 
vB u oc v 4 ). As the cloud becomes more and more optically 
thick, the lines become undetectable against the bright con¬ 
tinuum. 

Importantly for us, there can be an intermediate regime 
where the optical depth to Ha is large while the optical 
depth to He II 4686 is not so large, and so the Ha emis¬ 
sion can be somewhat sup pressed relative to t he He II 4686 
emission (see also lGaskell fc Roias Lobosll2014f ). This is the 
regime where our calculations generally find the values of 
£/ 4686 /£/Ha to be largest, and also where our earlier esti¬ 
mates based on the continuum suggest PSl-lOjh was at the 
time of the first spectrum (see f l A3I) . We emphasize that un¬ 
fortunately, these high fluxes and optical depths also lead to 
convergence issues for Cloudy, making Cloudy an ineffective 
tool for exploring this key region of parameter space. 

Now we compare these results for T = 2 x 10 s K with 
the results for other incident SEDs. Looking ahead to the 
other thermal SEDs shown in Figure IA2I we see that the 
structure is very similar, with slightly different temperatures 
leading to slightly different optical depths and line ratios. 
The maximum line ratios f/ 4686 /Lira are only 1.7 and 1.8 
for 5.5 x 10 4 K and 3 x 10 4 K respectively—but much of the 
region at large 4> and large n is unconverged, so perhaps we 
are missing larger line ratios there. 

Back to Figure EH the structure seen in Panels (a) - 
(c) in is similar to that in Panels (d) - (f), but the line ra¬ 
tio in the high-flux region of the figure is close to 1 rather 
than 0.4. The reason is that the AGN incident SED con¬ 
tains a luminous X-ray tail; for the same incident number 
flux <f>, the hard AGN SED is more energetic than the ther¬ 
mal SED, and so the gas is much hotter (by a factor of 
~ 100 or more). Recombination coefficients depend on tem¬ 
perature: at low temperatures they scale roughly as a oc 
T ~ 0,5 , while at high temperatures they scale as a oc T~ 1,5 , 
with t he dividing line ar ound T'/Z 2 ~ 10 6 K (iFerland et al.l 
Il992l ; lYerner fe Ferlandl Il996l l . Recombination coefficients 
for h ydrogen-like ions are related by a(Z, T) as Za(l,T/Z 2 ) 
(e.g.. lQsterbrock|[l989l : IVerner fe Ferlandl[l996ll . and thus at 
high temperatures, 04686 is higher relative to oh a than at 
low temperatures. This in equation (IA2I) explains the line 
ratios in the top part of Panel (a) in Figure EH 

We also note that for the thermal SEDs at the highest 
values of $ and n, the neutral fraction of hydrogen is high 
enough that the gas becomes optically thick to Ha, and so 
the line ratio 7/4686/Alfa can rise above unity. By contrast, 
the AGN SED heats the gas so much that collisional ion¬ 
ization is important as well as photoionization, and so the 
neutral fraction of hydrogen remains low, the gas remains 


© 2015 RAS, MNRAS 0 Q 0 .mi 26 l 
























22 L. E. Strubbe and N. Murray 


AGN SED with hard X-ray tail 



T = 2 x 10 5 K 



7 8 9 10 11 12 13 14 


log 10 n (cm 3 ) 



Figure Al. Panel ( a): Line ratio L 4686 /^Ha for a range of incident fluxes and number densities calculated with Cloudy C13.01, similar 
to iKorista Sz Goadl (|20Q4h ’s Figure 1 calculated with Cloudy v90.04: an AGN incident spectral energy distribution (see text), a column 
density of TVh = 10 23 cm -2 , and cosmic abundances. White squares rep resent models that did not converge. The maximum calculated 
line ratio is 2.8, almost a factor of 2 below the lower limit observed bv iGezari et al.l (|2012l ) for PSl-lOjh. Panels (b) and (c): Optical 
depths in the lines Ho; (b) and He II 4686 (c), for the same Cloudy calculation as in (a). The highest line ratio L^qsq/L n a comes from the 
region where Ha optical depths are high but He II 4686 less so, so that Ha; is relatively suppressed. Panels (d) - (f): Same as Panels (a) 
- (c) but for a thermal incident SED with temperature T = 2 x 10 5 K. The maximum calculated line ratio is 3.4. Qualitative similarities 
and differences between results for different SEDs are described in the text. 


optically thin to Ha, and the line ratio L 4686 /Xha remains 
close to 1. 

Because of the qualitatively different behaviors seen in 
the top row relative to the b ottom row of plots i n Figure 
[All we disagree with lGaskell fe Rojas Lobosl (|2014l Vs state¬ 
ment that results do not depend on ionization parameter 
(U = $/nc) or SED shape: Rather, it is crucial to perform 
calculations using the appropriate incident SED in order to 
properly capture how the temperature of the gas affects opti¬ 
cal depths and line emission. No X-ray observations of PSl- 
lOjh were taken around the time of the spectrum, so we can¬ 
not know the true shape of the incident SED that produced 
it; however, the late-time Chandra observation showed no X- 
rays. Further investigations of photoionization in tidal dis¬ 
ruption events must go beyond using only AGN SEDs. 

Finally, we point out that workers in the field occasion¬ 
ally offer the explanation that Ha is suppressed relative to 
He II 4686 because the hydrogen is “over-ionized.” We em¬ 
phasize that Ha emission is produced following ionized hy¬ 


drogen H + recombining with an electron (and He II 4686 is 
produced by He ++ recombining with an electron): a highly 
ionized gas cloud produces more Ha emission than a par¬ 
tially ionized one, not less. In optically thin systems, the 
number density of neutral hydrogen is irrelevant to deter¬ 
mining the luminosity of Ha emission (see H A2f) ; H° becomes 
important in this story only for absorbing Ha photons and 
suppressing Ha emission. 


A5 Problems with applying Cloudy calculations 
to PSl-lOjh 

Having understood the behavior of the line ratio L 4686 / Lu a 
in these suites of Cloudy calculations, we are positioned to 
ask how applicable these results are to PSl-lOjh, and what 
they can tell us about the nature of the event. 

A serious problem in applying these results to PSl- 
lOjh is that the Cloudy calculations assume that the gas 
is stationary, while the observed linewidth of He II 4686 
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T = 5.5 x 10 4 K 
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Figure A2. Same as Figure lAll but for thermal incident SEDs having temperatures T = 5.5 X 10 4 K (Panels (a) - (c)) and 3 X 10 4 K 
(Panels (d) - (f)), in line with the SED observed for PSl-lOjh. The maximum calculated line ratios l-.-icse/T hck are 1.7 and 1.8, for 
T = 5.5 X 10 4 K and 3 X 10 4 K respectively, well below the observed lower limit for PSl-lOjh. However, cooler temperatures lead to 
higher optical depths and more convergence issues for Cloudy, including the region of high flux and high density where largest values of 
^4686 / 7 h tv are most likely. 


~ 9000 km s -1 assures us that it is not. The lines in Cloudy 
are only thermally broadened, by ~ 10kins -1 for typical 
gas temperatures of ~ few x 10 4 — 10 5 K (appropriate for 
the thermal incident SEDs), or up to a few xlOOkms -1 for 
the most extreme cases with the hard AGN incident SED. 
This means that the line broadening should typically be 
~ 10 2 — 10 3 times greater than in the calculations. Cru¬ 
cially, this should diminish the optical depths in Ha and 
He II 4686 by this same factor of 10 2 — 10 3 : Doppler shift¬ 
ing due to large velocity gradients in the gas allows line 
photons to escape that would otherwise be trapped. We 
thus expect the gas to be close to optically thin to Ha 
(and He II 4686) even in the high flux, high density cases 
where our (stationary) Cloudy calculations showed Ha to 
be optically thick (see Panels (b) and (e) in Figures [All and 
IA2I) — and where we found the highest valu es of 1 /4686 /T hq- 
This should alter the radiative transfer llKorista &; Goadl 
l200d) . and likely reduce I/ 4686 /Ehq towards its optically 
thin value (~ 0.4 or ~ 1, depending on the gas temper¬ 
ature), even further away from the observed lower limit 


for PSl-lOjlS We believe this poin t casts the most seri - 


ous doubt on 


Guillodion ^fanukim fe Ramirez-Ruid (l2014h 


and iGaskell fc Rojas Lobosl ( 20141 4’s claim that PSl-lOih 
was produced from gas of cosmic abundances and could thus 
have been the disruption of a main sequence star. 

Another concern is that Cloudy requires conditions that 
are optically thin to Thomson scattering, r es <j 1. As ar¬ 
gued in H2.ll the observed blackbody continuum indicates 
that the gas giving rise to the emission is optically thick, 
and Thomson scattering is the dominant opacity for the ap¬ 
propriate temperature and density regime. We perform our 
Cloudy calculations by dividing the gas into two pieces: the 
Thomson-thin outer piece that we simulate with Cloudy, 
and the Thomson-thick inner piece whose continuum emis¬ 
sion (which we simulate simply by choosing various incident 
SEDs) irradiates the Thomson-thin outer piece. To investi- 


28 An additional concern is that broadening lines makes them 
appear fainter relative to the continuum; to agree with the obser¬ 
vations, photoionization calculations must result in a He II 4686 
line that is broad and nevertheless significantly brighter than the 
continuum. 
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gate spectra of TDEs properly, one should simulate the full 
extent of the gas all at once, including physics describing the 
luminosity generation (presumably) deep inside, and Thom¬ 
son scattering. Unfortunately, Cloudy is not an appropriate 
tool for such studies. 


A6 Conclusions and directions for future work 


iGezari et al. Hua measured a lower limit of L 4686 /Lhc* > 
5 in the spectrum of PSl-lOjh, taken 22 days before peak 
luminosity. In this appendix, we have presented Cloudy 
calculations aimed to answer the question of whether this 
observed line ratio could have been produced by gas of 
cosmic abundance; if not, PSl-lOjh presumably could not 
have been t he disruption of a solar-type star, as ha s been 
clai med by Guillochon, Manuki an fe Ramirez-Ruid d2014l l 
and lGaskell fc Roias Lobosl ( 2014fh 

We used Cloudy to calculate line ratios for a suite of gas 
densities and incident ionizing fluxes, for a hard X-ray-bright 
SED and three thermal SEDs, and cosmic abundances. From 
these, we draw the following conclusions: 


(i) Results can be qualitatively different (see upper and 
lower rows of plots in Figure lATl) depending on the shape 
of the incident SED, as X-rays heat the gas to high tem¬ 
peratures and reduce the ratio of He ++ recombination rel¬ 
ative to H + recombination. Previous work has used hard 
AGN SEDs, while the observations indicate a thermal SED 
at Teff ~ 3 x 10 4 K — 5.5 x 10 4 K with no evidence for X-rays. 

(ii) The highest line ratios we calculate are L 4686 /L^ a = 
2.8 for the AGN SED and 3.4 for the thermal SED at 2 x 
10 5 K. T hese bo th lie below the observed lower limit of 5 
found bv lGezari et al.l (I2OI2I ) , an d close to or b elow the lower 
limit of 3.7±25% inferred bv lGaskell fc Roias Lobosl d2014h . 

(iii) The highest line ratios are found for models having 
large ionizing flux 4> and large density n\ the gas is highly 
ionized throughout, but the high density leads to a high 
recombination rate and so a relatively large neutral fraction 
of hydrogen, and so a large optical depth to Ha. This can 
suppress Ha emission relative to He II 4686. 

(iv) However, the Cloudy calculations are performed for 
stationary gas. Large velocity gradients in the gas, as indi¬ 
cated by the observed He II 4686 linewidth ~ 9000 kms^ 1 , 
reduce line optical depths dramatically, and could reduce 
the line ratio 1/4686/L/Ha from ~ 3 at maximum down to 
~ 0.4 — 1, well below the observed line ratio. 

(v) Cloudy has convergence issues as the number of “re¬ 
solved” hydrogen levels increases, especially at high fluxes 
and densities. We suggest tha t these issues like l y affec t the 
validity of Cloudy results in iKorista fe Goad! (I2004T ) and 
ICaskell fc Roias Lobosl (l2014l l as well. 


We thus assert that, although these calculations are not 
conclusive, it is unlikely that the spectrum of PSl-lOjh was 
produced by gas of cosmic abundance, and thus it is unlikely 
that PSl-lOjh was the disruption of a solar-type star. For 
these reasons, this paper investigates instead the possibility 
that PSl-lOjh was the disruption of a stripped helium core. 

For future theoretical studies of tidal disruption spectra, 
we suggest Monte Carlo radiative transfer studies, such as 
those currently being undertaken by N. Roth et al. (personal 
communication). Important questions to examine include: 


(i) How do large velocity gradients affect the spectrum? 

(ii) How does optically thick Thomson scattering affect 
the spectrum? 

(iii) If a predicted line ratio 1/4686/L/Ha is large but the gas 
is optically thick at many wavelengths, are the lines bright 
enough to measure against the bright continuum back¬ 
ground? Can theoretical calculations reproduce not only line 
ratios but also line intensities relative to the continuum? 

(iv) How do results vary for different abundance ratios of 
helium to hydrogen? How much must hydrogen be depleted 
to produce the observed spectrum of PSl-lOjh? 


APPENDIX B: HE II LINES, CONTINUED: 
IONIZATION EQUILIBRIUM AND OPTICAL 
DEPTH 

In this section, we justify analytically several assumptions 
about He II lines made in fj4] 


B1 Helium ionization and equilibrium 


We solve the Saha equation to estimate the ionization state 
of helium. Assuming that the gas is purely helium, the ratio 
of doubly ionized to singly ionized helium is 


"He++ = J_ /2tt mJcT\ 3/2 e++ e —x He+ AT 

n He+ »el h 2 ) U He + 


(Bl) 


where [/ He ++ ~ 1 and [/ He + ~ 2 are the partition functions, 
and XHe+ = 54.4 eV is the ionization potential from He + 
to He ++ . In this section, our fiducial value for the num¬ 
ber density in the shell at the time of the first spectrum 
will be n S pect ~ 4 x 10 10 cm -3 (see eq. [5]). We find that 
for a wide range in densities of pure helium around this 
value, riHe ~ 10 8 cnU 3 — 10 12 cm -3 , the transition from al¬ 
most completely singly-ionized to almost completely doubly- 
ionized takes place steeply over a narrow range of tempera¬ 
ture T « (2 — 3) x 10 4 K. From photometric measurements, 
the tempe rature at the photosphere is inferred to be at least 
3 x 10 4 K llGezari et alJ[20l3) . Although the gas outside the 
photosphere at _R p h may not have exactly this temperature 
(since the gas is optically thin in the continuum), its tem¬ 
perature likely is not too different from this since the gas 
remains coupled to the radiation through line transitions. 
Therefore we expect that the gas is almost completely dou¬ 
bly ionized. For example, the fraction of singly ionized he¬ 
lium is n He +/n He ++ « 9 x 10 -3 at T = 3 x 10 4 K and 
n He ++ = 4 x 10 10 cm -3 , our estimates for the conditions 
in the shell at the time of the first spectrum. Furthermore, 
almost all of the singly ionized helium is in the ground state, 
since kT -C y H e+ • 

The time for the gas to reach this ionization equilibrium 
is the recombination time, 


(n e a He ++) 

^ n, 


10 


8 x 10 10 cm -3 


s, 


(B2) 

(B3) 


where a He ++ ~ 8x 10~ 13 cm 3 s ~ 4 is t he Case B He ++ recom¬ 
bination coefficient (IOst,erbrocklll989l . p.38) at T = 4x 10 4 K 
(which grows weakly with decreasing temperature). This 
timescale is much shorter than the expansion timescale of 
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the shell (~ 70 days to reach peak) or timescales of observa¬ 
tions (days). The gas therefore is easily in ionization equi¬ 
librium, and should be expected to produce He II emission 
line radiation. 


B2 Optical depth to He II photons 

Here we show that the region just outside _R p h is highly 
optically thick to helium Lyman a photons at 304 A (He + : 
n = 1 — > 2 ), and therefore capable of supporting a line- 
driven wind. We also show that this region is marginally 
optically thin to He II 4686 A photons (He + : n = 3 —»• 4), 
allowing us to see these photons as a broad emission line. 
We estimate optical depths in these lines by comparing with 
the (local) optical depth to Thomson scattering, dr es . 

The optical depth to AHeLa = 304 A photons along a 
differential path length is 

, n He+ n He+,l CHeLo , ,, 

dTHeLc = — — --- 1 - dr es , (B4) 

%e++ ^He+ 

where n He + 1 is the number density of n He + ions in the 
ground state, and a HeLa is the cross section to He La pho¬ 
tons. The thermally-broadened cross section to this transi¬ 
tion is 

_ 3^HeLo; ( 4m p ^ 1,/2 A IT)t . A 

HeL “’ th “ Stt U nkTj Alle ' 21 (B5) 

“ 2x10 "‘(3VToIk)’‘' 2 "' 1 “ 3x1 ° ,V 


the gas outside the electron scattering photosphere (where 
T es <j 1) is highly optically thick to He La, and therefore 
radiation pressure on this line transition can drive a wind. 

A similar calculation shows that the region is optically 
thin to He II 4686 A photons. The local optical depth to 
these photons is 

%e+ n He+,3 04686 j_ /D n\ 

074686 —-OT es , (H9) 

n He++ n He+ °T 

with velocity-broadened cross section 

3A4686 —1 A 7'T^1A^ 

04686,v = —-- V AlHe,43 (BIO) 

87T 

“ ? x ( 4500 kins - 1 )' 1jt ’ 

and AHe ,43 w 1 x 10 8 s _1 (see footnote 1291) . For a gas temper¬ 
ature T = 3 x 10 4 K, the fraction of He + ions in the n = 3 en¬ 
ergy level is n He + 3 /n He + = 3 2 e“B - 1 / 3 )x He +/kT w 7 x^ 0 -*^ 
far less than the fraction in the ground state. The optical 
depth to 4686 A photons is order unity just at the electron 
scattering photosphere _R p h, and falls outside, 


dr 46 86 


a ( n He+/%e++ \ 

( 9 x 10 - 3 ) 


{ n He+,3/ n He+ \ / Av 

( 7 X 10-® ) V 4500 km s~~ 


dr es .(Bll) 


That the gas is marginally optically thin to these photons 
outside the electron scattering photosphere justifies our ap¬ 
proximation of using the observed line luminosity to esti¬ 
mate the emission measure ( EM ) in 1 )4.21 


f Osterbrockj 1 1989 .. p.77) with the Einstein A coefficienlF^I 
AHe, 2 i ~ 1 x 10 lo s _1 . To calculate the full optical depth 
along a single line of sight, we must account for the fact that 
photons pass through gas at a variety of velocities whose 
spread (Av ~ 4500kins -1 , the observed half linewidth of 
the 4686 A emission line) is much larger than the thermal 
velocity (~ 20 kms -1 for T ~ 3 x 10 4 K). Calculating this 
in detail involves determining the velocity profile as a func¬ 
tion of radius, and integrating along different lines of sight 
through the expanding gas. However, as we don’t know the 
velocity profile well enough, we instead simply approximate 
that absorption is spread out evenly between wavelengths 
AHeLa(l —A v/c) and AHeLa, so the velocity-broadened cross 


B3 Line-driven wind mass-loss rate 


The rate of mass loss due to the line-driven wind M w i n a can 
be estimated by setting the rate of momentum imparted by 
photons to the gas to the rate of momentum carried by the 
outflowing gas. Since the majority of He + ions are in the 
ground state (ffiD, photons need to have an energy greater 
than or equal to the He Lya transition in order to exert 
pressure on the helium ions. We estimate the mass-loss rate 


A/wind 


L(hu > 40.8 eV) 

CVoo 


(B12) 


section is 

^HeLa,v 


O' HeLo:,th 


icr 


(s) 

V 


( _-_V 

\4500 kms- 1 / 


0.02 


(B 6 ) 

,4500 kms- 1 ; <JT ' ^ 

From m the singly ionized fraction is n He + /n He ++ ~ 9 x 
10 ~ 3 for T = 3 x 10 4 K and n He ++ = 4 x 10 10 cm -3 , and 
essentially all He + ions are in the ground state; the optical 
depth to He La photons is thus 

dTM , L , ~ 10 - (I f T ) dr. : (B 8 ) 
V 9 x 10 “ 3 / V4500kms- 1 / 


29 Einstein A coefficients come from “Persistent Lines 
of Singly Ionized Helium (He II),” National Institute 
of Standards and Technology Physical Meas. Laboratory, 
http://physics.nist.gov/PhysRefData/Handbook/Tables/heliumtable4, 
When we use these values in equations QUO and fBKfl) , we ac¬ 
count for the different fractions of ions in the different spin states 
according to their degeneracies. 


£>40.8 


2 x 10 43 erg/ s J \ 4500 km/ s 


Mq yr 


where we have normalized the velocity to that of the ob¬ 
served half linewidth of the 4686 A line, and L[hv > 40.8 eV) 
to a value corresponding to a radiation temperature of 
Teff ~ 5.5 X 10 4 K with R p h ~ 3 x 10 14 cm, appropriate for 
that temperature at the time of the first spectrum. However, 
because 40.8 eV is on the exponential tail of the blackbody, 
this estimate is very sensitive to the assumed shape of the 
incident SED. The predicted mass-loss rate would be 3 or¬ 
ders of magnitude lower for a blackbody with T ~ 3 x 10 4 K 
at the same time of observation, but could be much larger 
if the event produced significant X-rays (which is possible 
since no X-ray observations were made until almost a year 
after the optical discovery). 

Building on this very uncertain estimate, the wind’s 
^llWisity would be 


Pwind 


Af w i n d 

47rr 2 Uwind 


(B13) 
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a * 


gem 


/ Mwind \ 

f Rph N 

\ f ^wind \ 

l 0.02M© y r - 1 ) 

\3 x 10 14 cm/ 

' V 4500 km s- 1 / 


This density is ~ 10 2 times smaller than the density we 
inferred for the shell. However, we have neglected the effects 
of metal ions, including the lines that are prominent in O- 
star winds (C IV, Si IV, N V, etc., which are lower energy 
than He La), so this estimate should be considered a firm 
lower limit to the wind density. Spectroscopic observations 
of ultraviolet resonance lines in absorption (such as helium 
Lyman a, and the UV resonance lines just mentioned) would 
allow better estimates of the optical depth and mass in the 
wind. 


APPENDIX C: THE PHOTOSPHERE IN A 
SCATTERING-DOMINATED FLOW 

Our analysis in S I2.1I assumes that the emission from PS1- 
lOjh is blackbody, as do analogous calculations for super¬ 
nova observations (e.g.. lQuimbv et al.ll201ll : IChornock et al.l 
2014). In reality, the dominant opacity in the flow is Thom¬ 
son scattering rather than true absorption, meaning that 
the emission in the outermost layer of the flow is “modified 
blackbody” rather than true blackbody. Our preliminary 
estimates (using the photoionization code Cloudy) suggest 
that at peak emission, the absorption opacity is 10 -3 — 10~ 2 
times the scattering opacity, and larger earlier on (when the 
density is higher). The spectral energy distribution of a mod¬ 
ified blackbody is broader than that of a true blackbody, 
so the temperatures inferred from the observations may be 
(mild) underestimates. 

These effects reduce the assumed emission by 
(Kpabs/Kpscatp 2 (where K„,abs and K„, SC at are the absorp¬ 
tion and scattering opacities), though this reduction is 
slightly counterbalanced by the higher true temperature. 
This implies that the radius of the photosphere may ac¬ 
tually be a factor of a few larger than we estimated in 42.11 
and Figure [lj implying a larger velocity, larger mass of the 
shell, and larger kinetic energy. 

Specifically, for a photosphere that is larger than our 
previous estimate by / ~ 3, the velocity would be Vo ~ 
3000kms _1 , the mass of the shell would be M s h e ii ~ 
0.09 Mq, and the kinetic energy would be Ek ~ 8 X 10 47 erg. 
The kinetic energy would still be far below the energy of 
emitted radiation. The shell would carry ~ 60% rather than 
7% of the (originally bound) half of the mass of the disrupted 
core. 

We defer to future work more detailed investigation 
of the photospheric radius accounting for deviations from 
blackbody. 
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